“Calhoun 


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1998-03 


Single-frequency measurements using 
undersampling methods 


Chia, Eng S. 


Monterey, California. Naval Postgraduate School 
http://ndl.handle.net/10945/8971 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 
| (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist sia Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

ia) LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 

































































































































































































ah cer bheeage 
” mypses aft ra) 
moe sia ep owe bane ge 
siteene & t Ont ght D Pe oe Pd 
Geo We oe Se 
pan Spans Rnd 00 ater Jy BA 
eaten re palel ve 
seoptes pee ry te Ad Cheb ad ght 
' whoa gage it 8 repeat ge c eabpeeqied’ Tet thee 714015 82 25 
S re GAL GU AD Delt! . wo Pee LE a rina AGP asec 9 sk 
| yee pate te pokee f ws ou PO Saeed 1 ¢)% ro ed Par TL Cala 
Be ee bengal’ SAD Perey i r ote fad cys Conee Me hod a8 Pane 
Paige tote 8 ' Pa Wea ag ae ME 
® ‘ rg teehee Pemge te iP at tet PD eth ght ent 
' 1¢f ‘her ada ated end Shot wah mide tle 
Pear ir oy ee) Ll ee! 3 SPs stant Ad; 
“y Famers se cll uC olbge ge aime : mh) rst to G, Wap oth, Awd 
ie sy im mpepte aT geebencd Oey} 1d te Pgh 18 288 
a! J aetnety Be ArT Ta cid ee 9 fe Leda de LE St Aes te 
abe , rata Be ad OPA RIDE? 8 Ot F ' whe Op 
pec gspre ue een? rae Ue SSID ED ES oh Oe pissy 
: ve at, ay? BEV Pad Loe 2 pag ot ol a etl FD 
n af CONE pear 1 at ieee 
mi * 
PO finer gic 


bg tate wes 
er a i re reeey 













! ooyee™ tuee “ft 


py ieee 
, Uylig Gediayle 


Teoh, Ud aaalat Cir 
ther 
8h Gh itd olf pina oP 2 Ban as 
als oh? 


nitale@ed ICI? A angerd 
whe! 






>%, 
Bhd a Unf 




































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































nut wl 
gas Me 
a,amvett ¢ Cre ed ee 
& oe bythe git bowel al 
tibpene eed #4 r 
Ae ce arena Br i sce tpt we arwthed AQ? BF Os ah 
Amey iO rceqatetad Aedalads Te a oD Lome tat gh oP ts 2 A A 
‘ Mantes 8s gg eae Ving og wirath of tes AGL 8 at ' ad Pol Wale or Atal ad ‘Budisp rived Aen 
a wae +n dade a Ea opiate ko ME oft Ur7ny Leite ade dart! ad) ff Wivita dae 0. WHAT ebass 
pep wae adie Veg haroPa ge arms: AP Ne BNR e b0 ahi 1 MelreBald s Reet ad itie Be 
auda ah waited F 3 aaa wee he Pe hte ade Fahad vt ene ve cateat hit 9/1 PS"h 
ret Pa ttgt Fal & fi “3! ay oad ted af LS aap tat ‘ aie Se gee AT te Aptle tds solid 4 
am ’ 2 tip hearst ast gated aye gees a Aah Oya std ad G » i othe aR obra a) Vad 13a he Sek 0? 
cepa gete: at we at s peated tetes asee ye What ar re Tae ' pdadet at aAisar Wi Se 
“ Pr at Tg opr rey Ft Ue ad int OS 4 fran as" ere 
, ffs ie gh er both 4 yhey, Sod toh o Vane hp ates! hr ld oh Bp gb ae 
Yong ot : ahead iO oA at hat Oi gh Ateitad of b1tt 00 g8a! odie Wi apeteteDs 1A aot 
oe ow tip pra? etony 1h I at eda im Dae acd g Dab se lm? $! “ ree Pore er f 
: Has : Fan Fy) wg V8 FT Pap eta bylatet ad Abele Pod ne ab aeeg ers’ a8 pe oR VD BUG M ee ww sad at gore t ea Fe 
ae : Aue oe Tall Ate Giri at meat of Vapapatam aurat Haag ot Dee P ap ded cael s ve tee tt OA Gs Fad atl al og PMG a8 Nndhe Dard, Malek git As) 0 Pal oP FoF Bae MBS a Se 
' ' ' Cit iota ee a re Lh ee ater a os" ais epecests pabidgda ts Yen aa ah tot tate ik oe bt atalt: err, Ser heal ae 
Paks a theo ' Ser eag t Arora had Ed pdaten val at eet orale ga (env g ws! . erated G8 oF» ted ah aMtad tele’ 
. ea a o8 . ue 88 tn tate ene 10 1 Saber am res! Yel yeild 2 gree ene Hyer Nabot. Puyol Tabaladiinns eee 
. ' . om .tle ou tit # Peay ree oat, Ug ait in) Bw phere Cee ted ete td = te 
. Pe a EAR ee sa 14 a aah, press atl b¢— Say Ae tet ve ide parthed 90 tel vr, Mehae 
a ’ 8 .. ie wae ne oriole ee Pe TT bid ae ae ‘i Midnsads is oyeete remand 42 gael ZUeaS oe rad, 
' ' ine 1 ote ? get oC gtedagag eet Te VES) bee het he Atahedoe aFab Ma ae a Dk PVP AT Valente mA 
' . gant agll * ae degilaaes WN 1 @LbTA: HED bole ¢ th 2 detec ateawt 1b (Baas od Ab gt Sake aftal atte! eo 
. : = - F toon wt * a8p rppue sat perst we Sf a gd atlatig ‘att 614 oS Pee Sehaaye ghar v* 
F Pa F ' ‘ pal hee wera! 4 Prtetet ‘an te @ Lire"abs he) ot A oat gr oF Of 
oe ee aiate ae yrs | : OD Suh g ct Aint Daltiosdal ot How ah A wap ated: 
' ee ere eee ‘ peg PeNegcr tel ke it ap hatatetad: Votive of org ton ase tam nrg eiahet of ted oe eters 
, ' : 1 eq 8 wale me! atyt Usb ta tee Apiise 3 ory oe hw LeTlsteratan ewe searmteb ate teh as 1potet 
A a ' ‘ tote afte & GO Vy de bed Dro desiae svn) oie! NBed stp0 ae tre thd Pah 08 teint a Bete 
A ayy Ste aes re as CPi Tee r¥} ‘ aie . Gata ue woe gota #7 Tate oh Mess sate 
ae 1 rer leh anes va sf 4 abe ‘1 ‘ ty gts Owe : 1aede 49 Werk ID> BePrage: nee Face? statD cho Bet « 
; ‘ ' Veuns we oagead Wb aaeea ab atats “Dee rubs ab oa sanene 48: 84, Op ete tere SS ppebar Huma e 109 oo Bidar 
‘ i a ' Peer a tat Ete os atat wide lard ode: hal tnbab alt nat Ver tad as ot. ated Maths athe + eetate op habe ae Ashesto! 
6 Ee 2 atGee 9 vie fens Boel etal wi OT gee tafe amy ery! He Oapud- Dives ad Aare Ps ee Po eet 
a an Tee Ce PURAe ere | t ety tein Cee hada els ee ee te tL ada ae GF ah ab g hp A p Ba ahd ahn te 
; . . ' - ye teue 8 Of eee inarremecelsy art uh ahesfatast Brad ne ABS «agin s8 bi 0: a 
oe eee By tha: . ‘heres ulaesshint pape bs st Sar tye Fae 
" : oe 1 i : F me Lebcs tedales = Dales weve ete 
’ feelers: 98.5 ‘ ope ott ‘ 14 wine ‘ a beso oe Dog eared a? ard |More Oe 
: eee CRs PDE CY ee oboe a ete o UsSieiare abet ath Waly oo) th elVB Mead, me mm 0° abo” Ded oraraye® ac o% 
Pa o i oy ace . : 5 ie rousemntareasn ls " ign ’ CAC eahyees's ‘v Bact AND D 1a amulha (alread gua Bae 
ay * ee =e aie a 2g Svar (Soe wig 1 Set eaee VetatetgSrdeg Ve ’ fal pice: (ite tt S0,kult ets ee 
: ee s ' ter! Cee | 1 . Pr ee . a8 é PAUES Seeds ™ Aa awit: Perea ata dtr wal oops Th patad 
: . ' peeowt ' Pn ee .oe Nye Sue tele day of Ad av ite nse tetany? Be, J LOL) oh Petarat femal 
i ' . » ' tantao tt ‘ Dae LT phen Ce Le Te &e at Ngote Rep SOA OG atte Sedu arhoie he er te Cd 
‘ S AuteL Cio nieran® . en ne a te Setede PY oe ea vinete ta vt duptpectee see adh» Pc snd une 8 40h al nd De t= Nahe bs 
a ° . 4 : wim twi ae S ; Ro dpe bye et Tighten g de We teue Ay 4 © ge aE ye hehe DTP Oo 
ie moe . See ere pare oe fee B tuly ¢ Byoedtee : ped Hadvd ade AS gtehs tatate 
‘ eo: oe toe eee . Vipera Palode Tie 7o0 odode ee 
Are > Sine . ea ' oye Se oare : . ‘ e Ode a Gli debengtarade f shat hete’ svelotae’® . 
‘ ne, FR A i he a Teh od Pr ete re a Me gran Badyt ee rey, | dibe® ng 1.9 Bar? ~ Pe a ee 
7 1 Aun Woah . e ees Guyhe lores cote abetted at, ee te tiered ov hed 6) ol win Behe Dae 
‘a Per ra . Pa | ae Pant yee . ‘ . ane eh oe whats (4 Vy dg 8 A wbasda mata? © Kigte 
beaks o “8 ee aka) fan Sav Oem, gan " stgd st . te dacheh ote APrarh yas * 
; . ‘ Part . ‘ ae ‘ ve noise H : a @ a dodels Fea, ertotephy 
' ' ma eon ieege BS Pc ha Ae al oe ae pee et: omar} ws 
, es Se ae oat ; $ 3 bev adeb hese ta hee a seretaten ate ies, tee 
er a: oo agteaels Sate: oe Una ae noe . feseteras O18 at gta, Date eens af . ate sharin fosgne at 
‘ ’ . " « oe ee Am | 2 1 od devaft be wettqes © ws oP atehe art oh.- 8 afte O sm ba eta a) e aah bite ate 6 te be M 
Z . ‘ Teds pr @. ay aean?é ete Hat Maga Cate titers doh tye afotre Og? @ petted af adem Prt pee Cae ke Dalene Ae tn" sane at 
so ’ Pp . si ces vale . aCoteiites (s nye deta tas, 2 HC oe 8 Ae Pe aT he hu ‘p qttee bade sp) Fat" 
A 7 ou ' “a Pie ot 8 On ou «48 oe dtane a ete es ke) 
by ‘ ca e ef re at agree Ratt i ec oe 1 ‘ 1 det epee rtst o oats of 
eee : : ‘ ' ran Se Sere vei 8 SET kage ole dade ater eto ahalet Aveudgtetarat 
, x Fy . F ' og WG eaghul neil s, ‘ i go dag ebototg of Nedget oe? ateVatarsin ¢ 
a : 1 8 om ft (gee Carper ete. m3 1 Dg tates Db Seta tateds + nydal tqeeee be io Baad AA Ke eet ates Dan ay 
a ; ' : ' ‘: ; Be ; 7 r eae ’ a 1, %0) % by) Byte anh Ly estat pene te . 
. ‘ aed ore Hee ' ae stare’ Gem Oty holm: jaw ey be 
By ee Beye ie es , e o.8 ite vo Fae wg talte Bo : "ae rope eee 
; e . . ‘ es . Ya arr ee rey ie ae nc’ Retate% mm het regen whe le tale 
, ' : Perot wi . weneteoe ™e™® Raves saatete- 
aeosntt ' . hae eT ha 
vee ‘ vs Beer F ns whe » Myla wnle’ 
. : : : ‘ . : Peery 1 oon « pies : ie hate sapere 
7 . i oe ‘ o. plete ‘ aratigh eat yh . vaerhe’nsla Vater ate 1tarhe 
s « ', A ' os , ‘ Peery an. eyehia Jaap ta sue Mobehure Pate pation. tet tg fig haces ot Ve \s 
: i e oe bre ones ve Peo \ee Vorad ee hetaten a oles Tere eae caibe® my he 56, “oe 
m4 ' er hoe te one . tbe nage shelhal ctted bab ote webytatess 
7 1 ae A eo 8 ee © ‘ ‘ a A: 
° . ' . ‘ ' oh) gp le P ? ae 4 
. Ai ‘ . . an Se AU Mee ; 
. far oat ' ; a : r) ns te i ant og le Puedes ofl od Tah; entore! 
7 Pr a ? o e 1? re Serie ae ORPa ": 
. 1 rae) Pa . oo? Ln eT a fm eat! def 9 vial Wi i: ’ . we aha etalon 
« oy ‘ : . ’ a es t verges © tebhe'a? ST he ponies “ + 
‘ 2 . t ‘ foo * . ot 8 aut fedataredel Nee a tek ebay WE Ered at 
: i tos ye fab aeal at ne? obs @ he pontgen  Me%8 Tet 
« 1 . ” of cg" . ee %e tet & be fe 2% Ogata otatstel stete tar " ote otedatale 
. 1 . ae . s ! 4 otpeat 1 ' wifte Ye ated! 
' x ' Ber . ‘ oo 8 pees seace™ of e wetter 9 ole 
. et Pa 8 gtttg efar % Be oe 
". = : . ‘ ’ a 1 = ane eye ‘ 4bream ota aqt ole 
' 7 . ‘ is Oe ape ne ake stefet tote an de y weart 
. . . . . ‘ Je utye 0 tats 2 4 » > geet : ofa 
oF Ty : "> rhe gts o o% om Se Osa aee 
' a cr ee e= ° Fi car etucl Galea rei arity 9 eS eres 
. > . oA ae aia te ' ry ea i prs a 
2 a ." ‘ ° ests . i! at uietarargte's” 
-! . « * : “t . ‘ Ls! ad ath" « 
at . v4 a oa «8 
eo. . aC . ~*? @ . 1 ' Sage tte 4 
oe ' v« . . te « af % AG 
a ' ’ we i ot ” Sten aens Fi ; a ee al 
' : . 1 ' te 4 
a 8 ‘ . . « ’ « . ‘8 ' ¢ 
. ‘ « : 1 ig, ou? ‘ S..° vOprvtets FA Sg 
, ' . ° . . au ‘ah Y oe etude tard’ st 
r) r r . eo = ’ ot s VaqrrmabaGlats © 6 
. ‘ . ' e . 5 a a ee! “ sytt 
.* ° . 8 . aren ve Se r] ' . . om 
4 a . oe 8 Ly ‘ it pha ~ 
i] ‘ . e. a | es . = 
. e ° roy 8 ‘ ve 1 1° mt 
' , ’ i 
. « "4 ’ . . ot) ¢ ' vow 
Ae oe s optee « * % ’ « 
. . we . . Ee . . ’ t. a be 
1. ’ . a ¢ . « a = ie 
M4 : ‘ . a. (8 “ae s 
* ‘ . “a 8 e *' e°® 28 ‘ 
.? ’ 2 ' 4 Chane ' 
. ‘ ar acon ° ; ’ ae 
% eer Te oe. " Ch a es : t 
x y ‘ Poa . FECA © aguit s 
; Beats P . en Ce Crs) 
' r i F 78 eon apy 
. ry sue a F 1 et, TC ’ ae 2 8 
t 1 ‘ ' Pee ee ’ ' ? Alar ii in ta ‘ et 
. . ; ; 7 is a 
« ¢ Ph ‘ 5? ¢ e¢ ® 
. .« ‘ . 8 ' eta an » © . , 
: . . . ’ ota . es ‘ s 2 
’ Het 
u ' Bie) ° 
° ' poate 1 7 e*pae ’ » ea 
i. ’ . . 8 . i. «? ° or ‘ ‘ ‘ ) 
. . ‘ +? e. e ait oe 1 
; « . ‘, . ’ aes | eee he 
< = . : ‘ ° eee er ‘ ae f% 
: 1 . ’ . 8 1° ' Pye . apt 
: d on Ar ' * ' 1 re 1 
i = tel ' ae ae : coy gent ? fa eink ‘ 
1 . . rT 4 of ee . . . . e. La ' ‘ x 
. *». . ' 1? ’ ’ 
‘ ‘ : 5 ae . . Me rh sa . : Cert 
| . 1 o od] » # ‘ ' » rer ‘ 
: oo. . 1 eta’ ee ‘ om O8 ' ’ ' 
, : 1 : . : nly 8 ' , Pe eer 1 ft. 
. » ‘ ' ge se . ee" ry se ' 
. ee ? = ' ‘ ge fe 1 
' 2 . sy, ya : ‘ ry ‘eo fa : 
= : on ; ¥ a » eee oe a b + { ' 1 A 
e . ' . . . ' ' ' i ase sa a 5 
. « . a | e¢ 1 Lf ot (22 
z ' ' . . ' 1 . eo. *, . 8 9,0 -te* . ihe 
s 7 ' ‘ e . * te e 86 pa hele ae 
' ; ’ 1 tae . . 4 ' a al 2 ‘ . 
F ‘ 1 ' " , 7s oot n 
Ss oo ee ' . e «a! . 
. . s ‘ a te ’ . * 
. ' . eee ' ' 
oe : o¢ . . “ ue? 
. Pit oy ta oe ‘ os sor . ea | 
‘ Py 8 acu. 40 © s 1 
: *? ' . . toe ae ° $ » 
' « . id . Pr a © 
. ’ * ’ 
. t « . « . ‘ ry ve 
° . ’ . ab 
. * ‘ : ao se 
ef ® 8 . ot “i . 
. « ry . ‘ e,* , . 
: ana . i <a : 7 . Soe ae Beye 
ee = 
° . ‘ bod t 
. “« ‘ . 7 to a | ee F = 
oe Pa ‘ ‘ ; one + os 
. oat o's ' : qa? ' Mi ° 
. s . e . 4 » ° J , . ‘ 
° 1 
. Y . ve ; ang i i he 5 als % 
; F Par ° “6 s “ qbeteteSuate 28 
’ : ° Y : . Ng mang FES OFA 
[4 1 ’ «48 ’ . et ‘ Fart rah 
at . « : ‘ . or = yd a ; ai ", vytat 
7 7 ene 1: acne . me ye3°s 
. ' ' . ' ' a . ace e° a¢3 
. . 7 bat Vi caer mr of renee 
. r) 1 1 . oe 4 Py oes 4 4 +s 
e 1 oe ' ° 4 i] 4 fee aoe 
' ‘ 1? ! ‘ ’ a r F = Cris fc 
, or aah 3 se hr betetesg uy @arte%e 
eo fu aa o ¢* I ae t ° ra Tee! a 
‘ oe ' aa a’ : ae Zens teree OL” 
; : * ‘<r : ‘ toa ' U : fj egeee tae etten ake 
. . e : , x ‘ y ‘yea aeetas . ver aerate e 
‘ ' ‘ ' e -« : tee s ene ' eg AA RCE, I ee ete eens Pia! 
s ie S ees 7 ? : 4 +g ates oOre o%s eerars 
‘ " ae ' A ; Be ' A ‘. 8 ; etste toh f° } 
er or ‘ an) fase ' : : ; , SP ka see 
‘ ac <u + « » ¢ oe 8 ‘ ogee ’ ° ae soya! see 
A . ‘ ’ oe 8 Fi ' sebsey dee 
' of ' 8 ieee | d Stel sat asa eae 
' ‘ . ' a 4 . 8 1 . ss ¥ ¢: “as, . ALE Mataay, 
oa | ' 1 ‘ st ‘ ’ ' | eee bets aetna etutsa Vices 
‘ ' ‘ ia" ' . . ? ‘ wae Anas © “yen regret tetany 
» ¢ ' os , * 3 . ay —— sid eee Mees se are gyre eae 
: . , : , 7 on ae 1 6 “¢ Panes pea t-ie- te 
yy Ve ; ‘ ‘pos aa Ge : eee ” ee ba ay 
e 7 o 1 FI : ’ ; 1 ef " 4 Ea i ee ee « +4 Srnenetd 
; 1 ! i: a aha 3 ares 4a Ob" ys Og he rf 
ry ’ . . 1 1 eal Le ie eevee Ma Cocady' ate oe! 
' . ' 11 ' ’ ' . ' 1 4 ah ae syente oa made gas it ype: Ne wesrwste <* ap igt= 
on 1 ou - ’ Pia | ) ae ae fe i" wate * yn fn] ne gtsty's ® Pte Lea td Py st bd 
' . 2 © ms ? ' ' . ‘ " : . nae aie om” geet desvause't Pee eect eta? quaeregey ef pagrameee aot a“ f 
' ‘ 1 ' . A : ' 8 H oar ; . Pe tite merry oy eae 1 ra Tes gangte haps t wo nee Sa, ae sh ra a beeene’ 
U ' . P type Fi ear mm , ts . 1 7 i seqtete C Ot oneag hy ot «7A 4 *% 4 om? “ey ysgtar® 15) ey A fa 
= , ai te ' te Zaie , ! \ ; ooh ote teetes pio oy pen We dee gts"! Mo Ta oe ee hh Oe a 
i ‘ ‘ . - . opts tre ope at yey 8 ya Orethea a he | sume st Saat pcarerec & Hu 
‘ . a 8 ' Ld ' ts 0 9? ive we 4 ry +, ene tgtitt roan CPL hd be Feng g entered Lé 
. A a a : . ee ' 1 ©. : gre ath aus eu 2 aty ‘antee feel ‘a 4 4 o*t yee os) re de keene Ter 
: , ‘ : Dees ee he any a ae \ an ay er Lee tra te Need oe; ‘ Berg ty eery" 
: : : ’ 1 ' bao Y qgegte 14 : ia ; ~ 
. ’ ae me rtacnevapaveteiae ts “ a 1 4, rrr ete BN bt} pte Sy8 oe tad \ £8 oe lene 
: te on ' . ieee a " ‘ ef! 7 3.é eh kt etyte le SRO ig. v a npehetd ote 
: * ’ ' | cae | Paes fat i La . i eae . eft os S'Sty save orgie be ore yi Le nmirary ce eee’ 
as , Oe Ds age eee ee cece? daha ie weap etc a sSo ties pe starge din” ee Mery Pe Tereg te Coy” Tete sot epee HEE 
by etawe Svea ete a 5 srenelgial ve aide tal eeaeySrarttee at Paces gina? rig fee bebe 
#% 7 teu A ; Pa We oe 1 ' siti vns 1 antn be tT Ghp ge avec ately epiood yee. @ alee igiwaeee ryt at Py ie 2 ae 
: : ae _ : ot asere ean i ‘ : Pet a eee pro oe (eels g iyft gt ofy iy Ma Sot ' ee peo aged wren Me and te to Op Seiwa 19k grat Taek y+ It 1 ied See ee rgtereny eee? 
: '- 1 one Bab S Or a rett pea a ‘ yi Wait Saal LES Vtzie fo totnes! ie _ wth ne tg got Sf end Bene rete tt 9%. wees Sarr s se lereg saree at PGS 
: - : , 2 a) oe ‘ . rele A st pat orn PO ae'® yy ptyetsaeaeaty oth te arereZe 4° Se ads tgedesee of diss aer visa gtute” hk bat ota 2g Nee el 
: 7 ie . ‘ ney A Viangs ' mF mw vepeztgege Fe AeL HL NY ee 504 as bey fun tote iy ota ANk Sas piesa bet MP 
7 - va ' ‘ . a Trig ate sepe oe ‘ we p atgan at oRcae ¢ a csgditety $y nee tere sense ree vonprdnn a he4 
: , ‘ ‘ ‘ a “4 u Pd fe gh ante vee ot Path ] og pong ae eee po ba par ebritee 
' : : ae Cnt eer: nu sreh ov ghee Se be AURAL A 
‘ee 7 ; . 1 ta ' ' ’ y te toe t a ute San ne nay” erases 
' a . 7 ’ eee ee ol ae o fy Hh ta % ae ews big erty Pd eke ee aoe 
' 1 ‘ ‘‘ "4 oe @ te ! rpeues rom eT, Sar swe sh ir? Seiearey es we 
: ' ‘ ‘ ant Pe ae eae ae Sorate : Pr te ike Pica he oe oy ne cate see EAN 
e t > La 
Le ‘ eae en toes 6 fhed + i439 %8 : fe suey oF ee rte rene 
i ; ‘4 oben oe a eet | . ' «6 8% Pues que One AA 40 Pare sunyee teas e3 wong are a Ee 
, : : ‘ ' . se ie ye ermey Gea eg ed 9! CATE 
‘ ’ . ° eat f] ' » bo . . soe? ee be 
Pay a . “ ' . set : . ail A Sep Hetyes o th gtret one: 
. ‘ a ry ' . on 8 . : ataty ate H Gis hegre ra eae sdete 
: »e ; te eet! he ; ang te ZOE TES EL 
. s ‘ . a a é Ls itu ‘eve 4 oh 
. of LJ of 1. Pg hee" # 
¥ . . ge! od of y eds rte 
' 1 « - =. oe! * 
. wade ee a°ee 
aie ‘ : aie i. Seite Gyaate: sign pee 
‘ . tee “6° f ' ‘ i were a ik ah Sle aree 
' ‘ : eg ner rae oF ee ase i ete acarerd pa Wren ee nNOS 
--! ° E be = ee ae ee eee sige woop eit oe ree seasi te teeta ents WHR Et 
a - ae ’ = ¢ ‘ t ab ‘ A ore «'@ °R atge ete cule A pod hi Lae fat Re te TI Te ee ented veer 
: . ‘ 2 , epmetatg'y sant g 16m beg hs edt hp tt Segucseace: te m4 etd OED 
. . . yivle ates 4 > Saat &: furheard: Oe Wed ee ens 
: . 1 ' ' oe . ‘ Pa * s (a5 te. st t Ca ' LA eE IIL 1 ares 
' . ' is 2) Ns set : es rise ges tie Vebeneue oe vt 1 ga cgea tere wharyer geese 
. 1 eo . ‘ wee ae ra ou rege de? teouet pee 
: ' : ae ' J Re a 4 oye fe ree att tute everett Fe 
. : ae . ' F ia ToL aba. a cert aN 2 Cy ecere ree ed HONE r yet ee eas ree 
' ‘ #- . A A oe rg tao tye ote ace 36“ 4a en ey"8 OED 
‘ ' ' : ' : : : SPT ee ha ta tykshd Ne aca era hte Mr 
‘ v8 . = . » a ge PT | Fearne eed 8 ew 
peak : , e550 de Mpanresnirseere 2S 
. . ‘ ® . ‘, ey 
as? . et ‘ 1 
oe we 
Ache fqgr) 
". ay ive wld? 
ar rea a A 
he’ UA ies 
Wet ata 


DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY CA 93943-5101 








NAVAL POSTGRADUATE SCHOOL 
Monterey, California 





THESIS 


SINGLE-FREQUENCY MEASUREMENTS USING 
UNDERSAMPLING METHODS 


by 
Eng S. Chia 


March 1998 


Thesis Advisor: Phillip E. Pace 





Approved for public release; distribution is unlimited. 





| REPORT DOCUMENTATION PAGE —s| Farm Approvea 
| OMB No. 0704- 
| 0188 


| Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing ) 
instruction, searching existing data sources, gathering and maintaining the data needed, and completing and reviewing the collection of 
information. Send comments regarding this burden estimate or any other aspect of this collection of information, including suggestions for | 

| reducing this burden, to Washington headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson Davis 
Highway, Suite 1204, Ariington, VA 22202-43502, and to the Office of Managemen and Budget, Paperwork Reduction Project (0704-0188) 
Washington DC 20503. 





1. AGENCY USE ONLY (Leave blank) 2. REPORT DATE 3. REPORT TYPE AND DATES 


March 1998 COVERED , 
Masters Thesis 


4. TITLE AND SUBTITLE 5. FUNDING 
SINGLE-FREQUENCY MEASUREMENTS USING UNDERSAMPLING METHODS NUMBERS 
6. AUTHOR(S) 
Chia, Eng S. 
8. PERFORMING 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) ORGANIZATION 
Naval Postgraduate School REPORT NUMBER 








Monterey, CA 93943-5000 





9. SPONSORING / MONITORING AGENCY NAME(S) AND ADDRESS(ES) 10. SPONSORING / 
MONITORING 
Research and Development Office, Washington DC AGENCY REPORT 
NUMBER 


11. SUPPLEMENTARY NOTES 


The views expressed in this thesis are those of the author and do not reflect the official policy or position of 
the Department of Defense or the U.S. Government. 


MATLAB is a registered trademark of the MathWorks, Inc. 
12a. DISTRIBUTION / AVAILABILITY STATEMENT 12b. 


Approved for public release; distribution is unlimited. cele en 


13. ABSTRACT (maximum 200 words) = 


The objective of this study is to verify the Symmetrical Number System (SNS) | 
undersampling receiver architecture using software and to investigate implementation 
issues using digital signal processing (DSP) hardware. In the software design, a 
MATLAB program is written to determine a single sinusoidal input frequency using this 
receiver architecture. Each channel of the SNS undersampling receiver consists of a low 
speed ADC, a discrete Fourier transform followed by a constant threshold device to 
detect the signal’s frequency bin. The detected frequency bins are then recombined in a 
SNS-to-decimal algorithm to recover the frequency of the signal. Error rate performance 
in a Gaussian noise environment at the input stage is evaluated. In the hardware design, 
a sinusoidal waveform is digitized, discrete Fourier transformed and converted from the | 
SNS format to a decimal value using a single channel digital signal processor. 
Implementation difficulties and design issues are discussed. 


14. SUBJECT TERMS 45. NUMBER 
Symmetrical Number System, Symmetrical folding, Undersampling, Discrete Fourier OF PAGES 
Transform. 96 
16. PRICE 
CODE 
47. SECURITY CLASSIFICATION oF | 18: SECURITY CLASSIFICATIONOF | 49 security ctassiFi. | 2° 
THIS PAGE LIMITATION OF 
REPORT luntmecinea CATION OF ABSTRACT RESTRICT 
Unclassified nciassilie Unclassified 
UL 
NSN 7540-01-280-5500 Standard Form 298 (Rev. 2-89) 


Prescribed by ANSI Std. 239-18 





DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 


MONTEREY, CA 988}r8304a for public release; distribution is unlimited 


SINGLE-FREQUENCY MEASUREMENTS USING UNDERSAMPLING METHODS 


df 


Eng S. @hia 
Major, Republic of Singapore Airforce 
B.S., National University of Singapore, 1989 
Submitted in partial fulfillment of the 
requirements for the degree of 
MASTER OF SCIENCE IN ELECTRICAL ENGINEERING 


from the 


NAVAL POSTGRADUATE SCHOOL 
March 1998 


S Sm My, 
> 
Le 
14 
YC. 
NV 





DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
ABSTRACT 
Suse MONTEREY GA 93943-5101 


The objective of this study is to verify the Symmetrical 
Number System (SNS) undersampling receiver architecture using 
software and investigate implementation issues uSing Digital 
Signal Processing (DSP) hardware. In the software design, a 
MATLAB program iS written to determine a single sinusoidal 
input frequency using this receiver architecture. Each 
channel of the SNS undersampling receiver consists of a low 
speed ADC, a discrete Fourier transform followed by a 
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environment at the input stage is evaluated. In the hardware 
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transformed and converted from the SNS format to a decimal 
value uSing a single channel digital signal processor. 
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I. INTRODUCTION 


A. UNDERSAMPLING 


The digitization of a signal is usually governed by the 
Nyquist theorem where the sampling frequency is at least 
twice the signal bandwidth. The Nyquist theorem however, 
places a limitation only on the information that can be 
derived from a single set of digitized data [Ref. 1]. If the 
sampling frequency is less than twice the bandwidth of the 
Signal being digitized, aliasing and consequently ambiguities 
Occur. With additional information however, ambiguous 
frequency components due to undersampling may be resolved. 
Such information may come from, for example, trial sampling 
periods. Rader [Ref. 2] described how trial sampling periods 
can be used to recover periodic signals. The trial sampling 
period which yields the waveform of smallest variation is 
considered to be the correct period and the resulting 
waveform the correct waveform. 

Pace, Leino and Styer [Ref. 3] examined the relationship 
between the Discrete Fourier Transform (DFT) and _ the 


Symmetrical Number System (SNS) as a means of resolving 


Single frequency undersampling aliases. They.showed that the 
DFT encodes the frequency information of a signal in a format 
that is in the same form as the SNS. In addition, they 
proved analytically that aliases resulting from undersampling 
a single-frequency signal could be resolved using 2 or more 
channels. Bach channel in a SNS undersampling receiver 
contains a low speed ADC, a DFT and a threshold device to 
detect the input signal bin number in the frequency domain. 
The bin numbers from each channel are then recombined to 


resolve the signal’s frequency. 


Bu PRINCIPAL CONTRIBUTIONS 


First, this thesis verifies the SNS undersampling theory 
advanced by Pace, Leino and Styer [Ref. 3]. An algorithm is 
written and coded in MATLAB to prove the methodology and to 
show that the frequency of an undersampled signal can be 
accurately measured. The algorithm is also simulated in a 
Gaussian noise environment. Error rates for the different 
noise levels are obtained as a function of the signal to 
noise ratio. Since the Fast Fourier Transform (FFT) is not 
Suitable for computing DFTs in this application, alternative 


methods are suggested for real-time applications. 


Second, possible hardware implementation problems are 
investigated based on a Digital Signal Processing (DSP) 
plattorme Several problems were encountered: the need for 
stable sampling frequencies, large memories and alterna 
methods for computing DFT for’ fast response time. 
Integration into future EW receivers must take these factors 
into consideration. 

Undersampling offers several advantages [Ref. 4]. Dt 
allows the resolution of very high frequencies in EW 
receivers using low speed ADCs. This is especially so if 
several SNS channels are used. In particular, the use of 
undersampling in the design of receivers will reduce their 


cost and complexity. 


oe THESIS ORGANIZATION 


¢ 


In Chapter II, the relationship between the SNS and the 
digital frequency domain as wee by the DFT is examined as 
a means of resolving Single-frequency undersampling 
ambiguities. It shows how the frequency of a signal that is 
undersampled at two different sampling frequencies (two- 
channel) can be determined. In order to use lower sampling 


frequencies, the two-channel case can be extended to three or 


more channels. In particular the three-channel case is 
discussed. 

In Chapter III, algorithms for the two-channel and 
three-channel receivers are developed and coded in MATLAB to 
measure the frequency of an incoming signal. Each section of 
the software is explained in detail. Results are obtained 
based on different Gaussian noise levels. 

A feasibility study/design for the two-channel case is 
Carried out in Chapter IV using a DSP development kit. The 
Suitability of using a DSP platform and its associated 
problems are discussed. 

Chapter V states some conclusions and recommendations 


for future research. 


II. BACKGROUND INFORMATION 


A. INTRODUCTION 


Digitization of a signal is usually governed by the 


Nyquist criterion when the input Signal is bandlimited to 0 < 


f < £,/2 where f, is the sampling frequency. For higher 
frequencies (i.e. f > f£f,/2), the process of undersampling 
Gives rise to ambiguities. However, with additional 


information (or channels), the frequency components f > f,/2 
can be resolved. 

Pace, Ramamoorthy and Styer [Ref. 5] showed that the 
discrete Fourier transform (DFT) naturally encodes’ the 
frequency information of a Signal in the same format as the 
symmetrical number system (SNS). Consequently, aliases from 
undersampling can be resolved using this method. The theory 


é 


set forth is elaborated in [Ref. 3]. 


B. DISCRETE FOURIER TRANSFORM (DFT) 


Since all signals consist of sinusoids, for simplicity, 
a single frequency Sinusoidal waveform is used for analysis. 


Assume the Sinusoidal signal is 


x(t)=2cosa@t 


and after sampling 
x(n) =2 cos @ n. 


(2) 


The DFT of x(n) is given by [Ref. 6]: 


N-] 
X(Kepex(ew k=01,....N-1. 


n=0 


(3) 


Applying the DFT to x(n) results in a discrete spectrum where 
|X(k) |]? is the energy contained in the signal at each digital 
frequency @=27k/N. The spectrum X(k) has N indices with the 


digital frequency of each index given by: 


ozns,..2n 12 2p E*D ee 


Rs) oe!) 
N N N N N 


for N even 
(4) 


and 


fozn an SI og SOL an B—P pp OD for N odd. 


N Nie N N 
GS) 


The analog frequency corresponding to each index is obtained 
by multiplying each value by f,. Since signals with digital 
frequencies in the range % < ® < 2% are indistinguishable 
from signals with digital frequencies 0 < @ < 7, a digital 
frequency of each index can also be written as: 

] (N/2) On (N /2 -1) 3 


O2m—,...,27 : re for N even 
N N N N N 


(6) 


and 


oan .,..2nt NZ a5 LN/2] Ie <2 = for N odd. 
N N N N 


a) 
where Lx] is the floor function and represents the greatest 
integer less than or equal to x. Thus the spectrum X(k) 


resolves into N integer indices and incoming signals will map 


into unique bins: 


| ea | 
iS 
eo 

nN |Z 

N |Z 


me AW eee: ] for N even, 


sen a Ff for N odd. 
2 ape 


(8) 


(9) 


For example, for N = 5 (f, = 5 Hz and the sampling duration 
T, is 1 second), the output bins after the DFT are [0 1 2 2 
1] for input frequencies of [0 1 2 3 4] Hz. These DFT bins 
are repeated for higher frequencies as illustrated in Figure 
1. In this figure the abscissa corresponds to the incoming 
frequency and the ordinate corresponds to the bin into which 


the signal is resolved. 





is A As 
0 1 2 3 4 5 6 7 8 9 10 
Input Frequency 


Figure 1: DFT bin mapping for input frequencies f= 
Onto .10 for N =95 %ie= 5 hz samelingerem ececenc 


Cc. THE SYMMETRICAL NUMBER SYSTEM (SNS) . 


The SNS is composed of a number of pairwise relatively 
prime (PRP) moduli. The integers within each SNS modulus 
however, are derived from a symmetrically folded waveform. 
The symmetrically folded waveform corresponding to each SNS 
PRP moduli (m,), has a folding period equal to the modulus. 
The integer values within each SNS modulus are derived from a 
mid-level quantization of the symmetrical folding waveform. 
The formal definition of a symmetrical residue is given 
below: 

Definition: For an integer h such that 0 <h<m 

<= min, been 
(10) 
If this function is extended periodically with period m, 


that is, 


Xny nm = Xp 
(eae) 


where n € {0,+1,+2,..} then x, is called a symmetrical 
residue of (h+nm) modulo m. For m even, let x be the row 


vector 


10 


x 
= ae 
> 
peek 
N|B 
N]B 
Ps 
N 
1 
ae | 


(12) 


3) 
where Lx] again represents the floor function resulting 
in the greatest integer less than or equal to x. These 
two vectors consist of the symmetrical remainder 


elements x,, 0 < hh < m. 


D. RELATIONSHIP BETWEEN DFT AND SNS 


From the above, it is obvious that the DFT maps real 
Signals naturally into the SNS. That is, in Section C, if we 
let the modulus m represent the sampling frequency multiplied 
by the sampling time (i.e., f£,T,), then equations (12) and 
(13) are in the same form as equations (8) and (9) where N= 
f.T,. Thus the SNS provides a convenient framework for 
undersamping signal analysis. 

Table 1 displays the input frequencies and the resulting 


DFT bins for sampling frequencies 5 Hz and 6 Hz respectively. 


deal 


Input Frequency | DFT Bins 
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Table 1: Input Frequency and Resulting DFT Bins for 2 
Channel Example. 

The frequencies are resolved as described in equations (12) 
shale awe cit a me By considering two or more channels, it is 
possible to unambiguously resolve the signal frequencies in 
the dynamic range determined by the SNS. One method is to 
devise a look-up table similar to that shown in Table 1. 
However this method is inefficient for high frequencies; 
large memories are required. An alternative method is 
described below: 

Suppose there are r channels and the incoming frequency 
is within the dynamic range of the system. To carry out the 


SNS-to-decimal conversion, we need to solve £ = a,(mod m,) 


i be 


for i = 1,2,...r, where a, is the corresponding detected DFT 
bin for each @m,. The Chinese Remainder Theorem states that 


there is a unique solution modulo M=m,*m, ...*m,.. A standard 


method of solution is to find integers b, such that M*b,/m, 


1(mod m,) where 1 = 1,2,....r in which case the solution is f 


M*b,*a,/m, + M*b,*a,/m, + ... + M*b,*a,/m, (mod M). In Sections 
F and G below, examples are given to illustrate this 


calculation. 


E. DYNAMIC RANGE OF THE SNS 


Let m,,..., mM be xr pairwise relatively prime moduli, 
then the dynamic range, D (0:D-1) of a SNS system is given as 
follows: 

e “Tf all the modula” are odd, then the dynamic sang=sor 


the system is 


I=! =j+] 


De in| 7] m. r=)" m, 
2 hs 


(14) 


where j ranges from 1 to r-1 and m,,m,...m; range over 
all permutations of {1,2,3,..,r}. For example, for a 


two-channel case with m,= 5, me= 7, 


3 


cok ine 
D= min) Sa 
2 
or D=6. 
Poe nCee ehalleumcase with mm — Sate — 5) , —ee 


Wig 
D=—~min{m, sr IONIC, sien ela Ta =p m,m,} 


Ofte = 77. 


If one of the moduli (m,) is even, then the dynamic 


range of the system is 


D= ving m, +T] a, | 
l=2 


l=j+l 


(15) 
where j ranges from 1 to r-1 and m,,m,..m; range over 


all permutations of {2,3,...,r}. For example, for a 


two-channel case with m,= 6 m= 5, 


D= min] m. | 
Z 


Or D:= &. 
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For a three-channel case with m,= 8, m= 5, m= 7, 


. Mm, Mm, Mm, 
D=muin 5 eee apes es? 5 3 +m, 
or D=27. 


Clearly, the dynamic range of an SNS system with one 
even modulus iS superior to that using all odd moduli. 
Moreover, the greater the number of channels, the greater the 


dynamic range. 


F. THE TWO-CHANNEL CASE 

Figure 2 shows the block diagram of ae two-channel 
receiver architecture to determine a single frequency f. In 
this architecture the ADC sampling frequencies f., and f,, are 
relatively prime and T, = 1. The DFT outputs are thresholded 
to detect the frequency bins of the signal. The ———" 


frequency bins a, and a, are then used by the SNS-to-decimal 


algorithm to determine the frequency of the input signal. 


5 


frequency bins a, and a, are then used by the SNS-to-decimal 


algorithm to determine the frequency of the input signal. 











Anti-Aliasing 
Filter 


Figure 2: Block Diagram of a Two Channel Receiver 


Architecture. 
Let m, = f,, and m = f,, and suppose that the incoming 
frequency is within the dynamic range of the system. From 


Section D, we need to solve f = a,(mod m,) and f = a,(mod m). 
The two congruence equations, f = a,(mod m,) and f = a,(mod 


m) are solvable only if the greatest common divisor of m, 


and m, divides (a, - a,), a generalization of the Chinese 
Remainder Theorem [Ref. 7]. To solve for f, the diophantine 
equation 

p*m,+q*m, = (a, - a,) (16) 


16 


must be solved for p and f— is then calculated from the 
equation 

£ = a,+p*ire | (17) 
The code for this algorithm is shown in Appendix A. 

For example, for sampling frequencies 5 and 6, m, and m, 
have values of 5 and 6 respectively (T, = 1). If the signal 
is resolved into bins a, (= 2) and a, (= 1) after the DFT, p 
is found to have a value of 1 and g is found have a value of 


This 


| 
~] 


-1. Thus, the input frequency from (17) is 2+1*5 


can also be verified as shown in Table 1. 


G. THE THREE-CHANNEL CASE 


Figure 3 shows the block diagram of a three-channel 
receiver architecture to determine a single frequency f. 
Similar to the two-channel case, the ADC sampling treaenees 
£.., £,., and f,, are pairwise gS Tsp prime and T, = 1. The 
DFT outputs are thresholded to detect the frequency bins of 
the signal. The frequency bins a,, a, and a, are then used by 


the SNS-to-decimal algorithm to determine the frequency of 


the input signal. 


Ad 







Anti-Aliasinp 
Filter 


Figure 3: Block Diagram of a Three Channel 
Receiver Architecture. 


In the three-channel solution, let m, = f,,, m= f,, and 
m, = £,, and suppose that the incoming frequency is within the 
dynamic range of the system. We need to solve £ = a,(mod m,) 
and f = a,(mod m,) and £ = a,(mod m,). Using the Chinese 
Remainder Theorem and the Euclidean algorithm, the method of 


solution is to find integers b, such that M*b,/m, 


1 (mod m,) 
where i = 1,2,and 3 and M=m,* m * m. The solution is 
then f = + M*b,*a,/m, + M*b,*a,/m, + M*b,*a,/m, (mod M) where f 
is the frequency which falls within the dynamic range D of 


the system. 
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For example, let m, = 5, m = 6 andm = 7, so that M = 
210. and: =a. Suppose that the signal is resolved into 
bins a, (=@Ibe a, (= 2) anmd@aewe— 2) after the DFE. For the 


three-channel case the b, values must be found. Here, b, b, 


and b, are found to be -2, -1, and -3 respectively. Thus f 
+ 210(-2) (1) /5 = 210(-1) (2)/6 © 200(=3) (2)/7 mod(210) and we 
must choose the solution that falls within the SNS dynamic 
range D = 22 [10:21]. The correct combination £ = 84 -70 + 
180 mod(210) = 194 mod(210) . Although 194 is out of the 


dynamic range, 210 - 194 = 16 is in the dynamic range so that 


f = 16 is the correct frequency. 


FL 2 NOISE CONSIDERATIONS 


For a sinusoidal waveform, the Signal to Noise Ratio 


(SNR) is defined as 





(18) 
where P is the power of the signal and o’ is the noise power. 
Assuming a signal power of one, the noise power and amplitude 


are given by 
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nls 


(20) 


This o is multiplied by a normally distributed random number 
sequence of zero mean and unit variance and added to the 
input signal as noise. The simulation results are given in 


Ghapeer Til - 
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IIIT. SOFTWARE DESIGN AND RESULTS 


A. TWO-CHANNEL ALGORITHM 


The two-channel case was described in Chapter II. An 
algorithm was constructed based on Figure 2. The software 
given in Appendix A can be divided into the following 
sections: 

e Initialization. This section obtains all the 
parameters (number of iterations, input frequency, 
sampling frequencies, quantization levels) required. 

e Iteration loop. This section consists of a loop 
(with an initial count of zero) to count the number 
of errors. 

e Creation of Waveform. Based on the input frequency, 
a Sinusoidal waveform is created with noise added. ‘ 

e Sampling and Quantization. The waveform is’ then 
sampled at two different frequencies and quantized 
using a 14-bit ADC. 

e Windowing. A rectangular window operation of width N 
= f, * T, = £, (the total sampling/integration time is 


taken to be one)is carried out. 
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DFT Operation. A DFT is then carried out on each 
sample, taking only the first half of the DFT output. 
The formula used for the DFT process is a simple pair 
of nested loops. 

Bin Detection. A non-adaptive (constant) threshold 
bin detector is then used to find the bin with the 
maximum value for each DFT output. 

SNS-to-Decimal Algorithm. The SNS -to-decimal 
algorithm as described in Chapter II is then used to 


calculate the incoming frequency. 


A flow diagram of this algorithm is illustrated in 


Figure 4. The MATLAB code can be found in Appendix A. 
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Figure 4: Two Channel Algorithm 
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B. TESTING OF TWO-CHANNEL SYSTEM . 
To test the two-channel case (Sinusoidal signal without 
noise), the program is run with the following input and 


sampling frequencies shown in Table 2. 






Eee 





Dynamic Range 
[0:8 | tow input frequency 


100 97 0:145 Consecutive sampling 
A 
1040 547 1200 0:1146 Sampling frequencies 


meses! 12671) 12919 0:12794 High input frequency 





Table 2: Tested Input and Sampling Frequencies. 


For example, with input signal frequency at 7 Hz as 
shown in Figure 5, the sampled signals at 5 Hz and at 8 Hz 
are shown in Figures 6 and 7 respectively. The DFT output 
for the two samples are shown in Figures 8 and 9. The 
resultant bins of the first halves of Figure 8 and 9 are then 
supplied to the SNS-to-decimal algorithm to be converted to 


the input frequency of 7 Hz. 
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Progure 5. 





Input signal with frequency of 7 Hz. 
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Figure 6: Sampled signal at frequency 5 Hz. 
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Figure 8: DFT output with fsl=5 Hz. 
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It is found that if one of the sampling. frequencies was 
the same as the input frequency, the algorithm failed. This 
is because the resulting samples due to the same sampling 
frequency will consist of zeros. This problem can be solved 
by using at least two sets of sampling frequencies. Apart 
from this, the algorithm works well in this noise-free (high 


Signal-to-noise ratio) environment. 


C. SIMULATION PARAMETERS FOR TWO-CHANNEL CASE 


To obtain the error rates in a noisy environment, the 
two-channel software is run with the following parameters: 

e Number of iterations, num = 10000 

e Signal to Noise Ratio, SNRDB = -30 to 30 dB 

e ADC resolution, bit = 14 


e Input and sampling frequencies as shown in Table 3. 





Table 3: Input and Sampling Frequencies. 


D. RESULTS FOR TWO-CHANNEL CASE 


The results obtained are shown in Figure 10. 


S18. 


Error Rate % 


Figure 10: Error Rates vs. 


Srl: 


f=9 Hz 
f=90 Hz 
f=900 Hz 
f=9000 Hz 





SNR for two-channel system 


The following observations are made: : 


As expected, the error rates improve as the SNR 
increases. A tradeoff between SNR and error rate is 
required. 

Improvements in error rates were obtained when higher 

Frequencies were used. This is because at higher 

frequencies, higher sampling Frequencies are 

required. This leads to a higher N-point DFT (higher 
gain) which is less affected by noise. 

However, at higher frequencies, the time taken to 

compute the DFT was much longer. To reduce the time 

taken, the following methods can be implemented: 

e If N is highly composite (factorable into powers of 
many small prime factors, preferably primes < 10), 
use a "mixed-radix" FFT implementation. 

e If N is prime, or contains very large prime 
factors, use the "chirp-z" transform. 

e Use three or more channels in the receiver. A 
three-channel receiver has a higher dynamic range 
for the same magnitude of sampling frequencies. 


For example, a two-channel receiver with sampling 
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frequencies 6 and 7 has a dynamic range of [0:9] 
while a three-channel receiver with sampling 
frequencies of 5, 6 and 7 has a dynamic range of 


fOr |e 


THREE-CHANNEL ALGORITHM 


The three-channel algorithm is similar to the two- 
channel algorithm as shown in Figure 11. The MATLAB code 


can be found in Appendix A. 


a 






Enter num, f 
sl, f62 and fs3 





Error 


Count=0 


Create 


sinewave 


Add 


Gaussian Noise 





fsl f£s3 















Sample and add Sample and add 






Sample and add 


Quantization Quantization Quantization 









Noise Noise Noise 






Windowing 
and DFT 


Windowing 
and DFT 


Windowing 
and DFT 










Find Largest Find Largest Find Largest 






Bin Bin Bin 


SNS-to-Decimal 


Algorithm 


Increment 


Error Count 





Figure 11: Three Channel Algorithm 
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F. TESTING OF THREE CHANNEL ALGORITHM 
To test the three-channel case, the program is run with 
some of the following input and sampling frequencies in Table 


4. 


Remarks 


Low input frequency 


Consecutive 
sampling 
Frequencies 


Sampling 
Frequencies far 





Table 4: Tested Input and Sampling Frequencies. 


Apart from the anomaly discussed in the two-channel 


case, the algorithm works well in this noise-free (high 


signal-to-noise ratio) environment. 


G. SIMULATION PARAMETERS FOR THREE-CHANNEL CASE 
To obtain the error rates in a noisy environment, the 


three-channel software is run with the following parameters: 


e Number of iterations, num = 10000 


35 


e Signal to Noise Ratio, SNRDB = -30 to 30 dB 
® ADC reselmiezon, bit l=—14 


e Input and sampling frequencies as shown in Table 5. 





900 fat | ae 
Table 5: Input and Sampling frequencies. 


RESULTS FOR THREE-CHANNEL CASE 


The results obtained are shown in Figure 12. 
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Error Rate % 





Figure 12: 


Error Rates vs. 


f=9 Hz 
f=90 Hz 
f=900 Hz 
f=9000 Hz 


SNR(dB) 


SNR for three-channel system 


om 


« 


Comparing the two-channel and three-channel cases, the 


following observations can be made: 


The three-channel system is much faster than the two- 
channel system since the DFTs required are smaller 
due to the smaller sampling frequencies. 

However the results for the two-channel system with 
noise are better. For example to achieve a 
relatively error-free system for a frequency of 9000 
Hz, the two-channel case requires only -22 . dB. 
However, the three-channel case requires at least -4 


Gils) 


Sue) 


IV. HARDWARE DESIGN AND FINDINGS 


A. INTRODUCTION 


In the last chapter, the instantaneous measurement of 


frequency using the SNS-to-decimal algorithm was verified. 


There is a need to investigate the implementation of the 


algorithm in hardware. Digital Signal Processing (DSP) 


hardware was selected for the following reasons: 


A major part of the algorithm is the processing of 
DFTs which is a digital signal processing task well 
suited to be carried out by DSP hardware. 

DSP hardware provides a fast way to implement the 
algorithm. The DSP development kit is easy to learn, 
program and simulate. It is ideal for this 
application to investigate hardware problems. and 
limitations. 

Cost consideration: the development kit plus tools 
cost $1500; 


EW receivers are likely to incorporate DSP hardware. 
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S. TI TMS320C54X DSP DEVELOPMENT KIT “ 


The TMS32054C54x DSKplus [Ref. 8-12] is a low cost DSP 
Starter kit that gives a designer a working knowledge of DSP 
ae to build DSP based systems. The development kit 
contains a stand-alone application board that can be 
connected to the PC. It executes code in real time at 40 


MIPS while the Windows-based debugger analyzes it line-by- 


line, displaying internal DSP register information in 
multiple windows and in real time. It has an Analog 
Interface Circuit for the input of signals. The board’s 


communication interface enables the creation of C54x DSP code 
and host PC code. Moreover, the hardware enables the use of 
expansion slots for adding memory, peripherals such as 
interface logic, other DSPs etc. The developed code can 
eventually be loaded into a resident DSP processor, which may 
be part of a EW receiver architecture. Figure 13 shows a 
block diagram of the development kit. A more detailed 


description of the kit can be found in Appendix B. 
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ADC DSP DAC 


Parallel Portto PC To Other Interface Ports 


Figure 13: Block Diagram of DSP Hardware. 


ae SOFTWARE 


The software for the two-channel case described in 
Chapter II (Figure 2) is written using the DSP development 
el. The software (found in Appendix C) is coded in ‘C’ 
language/assembly language and converted to the C54x assembly 
language (if required) prior to execution: 

e Firstappl.c/Firstapp2.c. These two programs poll the 
input channel and sample the input signal at the two 
sampling frequencies respectively. 

e Hostappl1.cpp/Hostapp2.cpp. These two programs 
display the samples of the signals based on the two 


sampling frequencies and save the data in text files. 
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e Main.c. This program reads the data, executes the 
DFT, obtains the largest values for the two channels 
and then carries out a SNS-to-decimal conversion. 


These programs were run individually and consecutively. 


D. TESTING AND RESULTS 


Using data generated by MATLAB, the main program was 
tested successfully in the development kit. The programs 
were then run with an input frequency of 126 Hz and sampling 
frequencies, 125 Hz and 128 Hz. Results obtained were 
intermittent 1.e., correct results were not always obtained. 
A frequency counter and an oscilloscope were set up and it 
was found that the sampling frequencies were not stable. 
Testing with different frequencies did not improve the 


results. 


E. PROBLEMS 


Several problems were encountered elbhant ale: the 
investigation: 
e Stability of Sampling Frequencies. The development 


kit carries out frequency division of the master 
oscillator to obtain the sampling frequencies. 


Unfortunately, the crystal oscillator has a 
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resolution of 5 =JG33Hz. This iS unacceptable as a 
shift of 1 Hz in the sampling frequency will cause 
erroneous results. Moreover, the fact that the 
sampling frequencies are factors of the oscillator 
frequency and that they need to be pairwise 
relatively prime severely limits the choice of 
frequencies. A possible solution is to obtain the 


sampling frequencies directly from stable signal 


sources. 
Dia For higher frequencies, the execution of the 
DFT takes a long time. Several solutions were 


suggested and discussed in the previous chapter. 

Memories. Insufficient memory error messages were 
encountered when high frequencies were used. The 
same messages occurred when attempts were made to run 
the routines together. More memories and/or more 


efficient DFT algorithms are required. 


43 





V. CONCLUDING REMARKS 


The main contribution of this thesis is the verification 
of the relationship of the DFT to the SNS to resolve 
undersampling ambiguities and the investigation of hardware 
implementation issues uSing a DSP platform. Error rates for 
different SNR are also obtained. 

The use of undersampling technique using the SNS to 
measure frequency is a viable method to implement in a EW 
receiver architecture. However, the need for faster DFT 
computation and stable sampling frequencies must be taken 
into account before they can be considered for incorporation 
into EW receivers. There is also a trade-off between the 


number of channels and SNR. For faster response, a multi- 


channel case is recommended; but a higher SNR is required. - 
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APPENDIX A 


MATLAB CODE FOR SOFTWARE ALGORITHM 


ov? 


ove 


Thesis Project 


oe 


Two Channel Receiver 


oe oe oe 


oe 


clear all; 
S Initialization 


num=input ('Enter Number of iterations:'); 
fF=input ("Enter Input eeequency : '); 

fsl=input('Enter Sampling Frequency 1:'); 
fs2=input ('Enter Sampling Frequency 2:'); 


fp1=fopen('c:\matlab\bin\thesis\result.dat','at' 


%* Quantization levels 


bit=14; 
qnlevel=2“bit-1; 
g=2/qnlevel; 


for SNRDB=-30:2:30 


count=0; 


for i=l:num 
SNR=10* (SNRDB/10) ; 
sigmasg=1/2/SNR; 


t=(0:.001:1); 
Sig-Sin(2*putt ae); 
til—I7 fst: | tote 


noisel=sqrt (sigmasq) *randn(1, length ( 


ADCsigi=sin(2*pi*f*t1)+noisel; 
ADCsigl=fix (ADCsigi1/q) *q; 
t2=1/fs2:1/fs2:1; 


noise2=sqrt (sigmasq) *randn (1, length ( 


ADCsig2=sin(2*pi*f*t2) +noise2; 
ADCsig2=fix (ADCsig2/q) *q; 


%figure (1) 
sSubploe(s, 11) > pler(Ecna) 
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oe 


ov? 


oe? 


oe 


oo 


oe 


oP oi? 


oe 


oe 


is 


of of ol? 


of ct 


oP 


1) ane 


2) ) si 


Note: The sampling frequencies should be relatively prime 


ov? 


Number of iterations 
Frequency of signal 

Sampling frequency 1 
Sampling frequency 2 
; * Store results 


ae ol? 


~—— o\? 


No. of quantization levels 
quantization size 


Set Signal to Noise Ratio 
from -30 dB. to 30 db 


Error Courne 


Convert to non-dB units 
Noise normalization assuming 
signal power of 1 


signal 

first ADC 

noise 
digitized signal 
quantized signal 
second ADC 

noise 
digitized signal 
quantized signal 


stitle(*Figqure i. Plot of signal x 

sxlabel (‘Time’ ) 

sy label (‘Amplitude’ ) 

4Subplot (3,1, 2 oot (t1, ADCs) 

stitle('Figure 2. Plot of sampled signal (sampling frequency 1) 
plus noise') 

sxlabel ('time') 

ylabel ('magnitude' ) 


SUD ote, 3) ,plot (t2,ADCsig2) 

6title('Figure 3. Plot of sampled signal (sampling frequency 2) 
plus noise') 

sxlabel ('time') 

sylabel ('magnitude' ) 


%* Window operation 
% Assume rectangular window 


oe 


size of window is fsl 
size of window is fs2 
windowed sampled signal 1 
windowed sampled signal 2 


winsizel=fsl1; 
winsize2=fs2; 
winsigl=ADCsigi (1:winsizel) ; 
winsig2=ADCsig2(1:winsize2) ; 


oP o\? 


oe 


* DFT Operation 
DFTsigi=abs (fft (winsigi,winsizel) ); 
DFTsig2=abs (fft (winsig2,winsize2)); 
DFTsigla=DFTsigi (1:length(DFTsigl)/2 +1); % Taking half of image 
DFTsig2a=DFTsig2 (1:length(DFTsig2)/2 +1); * Taking half of image 


%figure (2) 

% Plot to locate position of maximum value 

% Note that due to MATLAB (which cannot have a zero index, the 
actual location is one less 

ssubplot (2,1,1), stem(DFTsigla) 

6title('Figure 1. DFT plot of signal with sampling frequency i') 
éxlabel('frequency bins') 

sylabel ('magnitude' ) 

subplot (2,1,2),stem(DFTsig2a) 

title('Figure 2. DFT plot of signal with sampling frequency 2') 
éxlabel ('frequency bins') 

sylabel ('magnitude' ) 


% bin detector 


oe 


yl, y2 are locations of max values 
Note that due to MATLAB, the 
actual location is one less. 


[i, yl] =max (DFTsigla) ; 
[j,y2] =max (DFTsig2a) ; 


oP 


oP 


al=yl-1; 
AZ —v ee, 
% SNS to Decimal Algorithm 


50 


end 
error 


o\e 


To solve for f==ai(mod mi) (where "==" indicates congruence and 
mi are pairwise relatively prime), the Chinese Remainder Theorem 
states that there is a unique solution modulo M=m1*m2....mr. 


oe 


o\e 


oe 


A standard method of solution is to find integers bi such that 
M*bi/mi==1(mod mi) where i=1,2,....r in which the solution is 
f==M*bl*al/ml + M*b2*a2/m2 + ... + M*br*ar/mr (mod M) 


oe 


oe 


o\e 


For a 2 channel case, i.e. i=1,2, 


% m2*b1l == 1(mod m1) 
% ml*b2 == 1(mod m2) 
% £f == al(mod m1) 


oe 


f == a2(mod m2) 
f == m2*bl*al + ml*b2*a2 (mod ml1*m2) 


o\? 


%* Given ml(sampling frequency 1) and m2 (sampling frequency 2), to 
* find bl and b2, the congruence equation is transformed to a 
* Giphantine equation and solved using the Euclidean algorithm: 


% m2 *bl may = 1 
% Mt *%o2 -nizsy2 = 1 
* The above two equations can be combined into 
% m2* bie mixb7—) 1 


* bl and b2 are solved by the function "lde.m" which is called by % 
"glde ; m" f 


f == al(mod ml) and f == a2(mod m2) is solvable only if the 
* greatest common divisor of ml and m2 divides (a2 - al). 

To solve for f, r from the diophantine equation 

6 r*ml+s*m2 = a2 - al must be solved. 

¢ r is obtained from "glde.m" and £f is calculated by the 

% equation f = al+r*ml | 


1diff=a2-al; 
r= Olde noiee ts lear. ) > 
freq=abs (al+r*fs1) ; 
% Count the number of correct results. 
if freq== 
COUNE=CounE+.; 


end 


= 1-count/num; 


%* Write results to file 


xl=fprintt (ipip'sd +d sd td +d td\n', £, fsl, £S2;/9SNRDB, Mum, cEeeE; 


plot (SNRDB,error, 'y+') 
title('Error Rate vs. Signal to Noise Ratio') 


al 


xlabel ('SNR(dB) ') e 
ylabel('Error Rate 7%') 

hold on 

end 

fclose(fpl) ; 


* To calculate the dynamic range 


if rem(tsl, 2) —— %* To check whether fsl is even 
DR=fsl1/2 + fs2; 

elseif rem(fs2,2)== 
DR=fs2/2 + fsl; 


oe 


To check whether fs2 is even 


else 


o\? 


DR=.5* (fsl+fs2); fsl and fs2 are odd numbers 


end 


o\? 


o\? 


Thesis Project 


oP ol? 


Three Channel Receiver 


oe 


o\? 


clear all; 
close 


t=} 


~ Initialization 


num=input ('Enter Number of iterations:') ; * Number of iterations 
f=input (‘Enter Input Frequency:'); * Frequency of signal 

fsl=input ('Enter Sampling Frequency 1:'); * Sampling frequency 1 
fs2=input ('Enter Sampling Frequency 2:'); * Sampling frequency 2 
fs3=input('Enter Sampling Frequency 3:'); % Sampling frequency 3 


fpl=fopen('c:\matlab\bin\thesis\result.dat','at'); 
* Store results in file for later processing if required 


* Quantization levels 


% bit=input('Enter ADC resolution:'); 
bit=14; 
qnlevel=2“bit-1; 
q=2/qnlevel; 


oP 


No. of quantization levels 
quantization size 


o\? 


oe 


for SNRBB=-30:2:30 Set Signal to Noise Ratio from -30 


aB to 30 dB 


oe 


count=0; 
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for i=l:num 
SNR=10~* (SNRDB/10) ; 
sigmasg=1/2/SNR; 


oe 


Convert CoO non-dB units 
Noise normalization assuming 
signal power of 1 


oe 


oe 


ov? 


first ADC 
noise 

digitized signal 
quantized signal 


tl=1/fs1:1/fs1:1; 
noisel=sqrt (sigmasq) *randn(1,length(tl1)); 
ADCsigl1=1000* (sin(2*pi*f£*t1)+noisel) ; 
ADCs igl=fix (ADCsigl1/q) *q; 


é 


oP ol? 


oo 


oe 


second ADC 
noise 

digitized signal 
quantized signal 


t2=1/fs2ats2:-1- 
noise2=sqrt (sigmasq) *randn(1,length(t2) ); 
ADCsig2=1000* (sin(2*pi*f£*t2)+noise2) ; 
ADCs ig2=fix (ADCsig2/q) *q; 


oP ol? 


oe 


oe? 


third ADC 
noise 
digitized signal 


t3=1/fs3:1/fs3:1; 
noise3=sqrt (sigmasq) *randn(1,length(t3)); 
ADCsig3=1000% (sin(2*pi*f*t3)+noise3) ; 
ADCsig3=f1ix (ADCsig3/q) *q; 


oe 


ov? 


figure (1) 


ssubplot (3,1,1) , pleat, ADCsigl (1:fsiee 
Stitle('Figure 1. Plot of sampled signal (sampling frequency 
1) plus noise') 
xlabel ('time') 
sy label ('magnitude' ) 
Beubplot (3,1, 2) , plot (t2, ADCsig2 (1-£s25 
stitle('Figure 2. Plot of sampled signal (sampling frequency 2) 
plus noise') 
6Xlabel ('time') 
tylabel ('magnitude' ) 
ssubplot (3,71,3) plot (t3 ,ADCsig3 (17£s3)) 
stitle('Figure 3. Plot of sampled signal (sampling frequency 3) 
plus noise') 
sxlabel ('time') 
sylabel ('magnitude' ) 


oe 


Window operation 
Assume rectangular window 


oP 


oo 


size of window is fsl 
size of window is fs2 
size of window is fs3 
windowed sampled signal l 
windowed sampled signal 2 
windowed sampled signal 3 


winsizel=fsl; 
winsize2=fs2; 
winsize3=fs3; 
winsigl=ADCsigl (1:winsizel) ; 
winsig2=ADCsig2 (1:winsize2) ; 
wins ig3=ADCsig3 (1:winsize3) ; 


oo 


ae ol? 


oP 


oe 


* DFT Operation 
DFTsigl=abs (fft (winsigl,winsizel) ); 


a5 


DFTsig2=abs (fft (winsig2,winsizez2) ); 
DFTsig3=abs (fft (winsig3,winsize3) ); 


DFTsigla=DFTsigl1(1:length(DFTsigl)/2 +1); 
DFTsig2a=DFTsig2 (1:length(DFTsig2) /2 +1); 
DFTsig3a=DFTsig3 (1:length(DFTsig3)/2 +1); 


oe 


Taking half the image 
Taking half the image 
Taking half the image 


o\e 


oe 


figure (2) 

Plot to locate position of maximum value 

Note that due to MATLAB (which cannot have a zero index, the 
actual location is one less 


subplot (3,1,1), stem(DFTsigila) 

Stitle('Figure 1. DFT plot of signal with sampling frequency 1') 
sxlabel ('frequency bins') 

sy label ('magnitude') 


subplot (3,1,2),stem(DFTsig2a) 

Stitle('Figure 2. DFT plot of signal with sampling frequency 2') 
xlabel('frequency bins') 

sylabel ('magnitude') 


subplot (3,1,3),stem(DFTsig3a) 

Stitle('Figure 3. DFT plot of signal with sampling frequency 3') 
%xlabel('frequency bins') 

sylabel ('magnitude') 


% bin detector 


oe 


[i,y1] =max (DFTsigla) ; yl, y2 and y3 are the locations of 
maximum values 
Note that due to MATLAB, the actual 


location is one less. 


of ole 


[j,y2] =max (DFTsig2a) ; 


oe 


{k,y3] =max (DFTsig3a) ; 


al=yl-1; 
a2Z=y2-1; 
a3=y3-1; 


* SNS to Decimal Algorithm 

bi=lde(fs2*fs3,fsl1); 

b2=lde(fs1*fs3,fs2) ; 

b3=lde(fs1*fs2,fs3); 

cl1l=b1*fs2*fs3; 

c2=b2*fs1*fs3; 

c3=b3*fsl1*fs2; 

EHO lola @l+d2*C2+a5*C3;al*clt+a2*c2-a3*C3; ai*cl-az*ce2+a3*c3 ; 


al*cl-a2*c2-a3*c3;-al*cl+a2*c2+a3*c3;-al*cl+a2*c2-a3*c3; 
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~al*cl-a2*c2+a3*cC3;-al*cl-a2*c2-a3*c3]; 
freqmat=rem(freqmat,fs1*fs2*fs3) ; 


for i=1:8 
if (freqmat (i) <0) 
freqmat (1) =freqmat (i1)+fs1*fs2*fs3; 
end 
end 
freq=min(abs (freqmat)); 


%* Count the number of correct results. 


1£ freq==r 
count=count+l1; 
end 
end 
Si0r = i-COuUnt/ num: 


% Write results to file 


xl=fprintf(fpl,'td td td sd td td td\n', £, £sl, £s2, £s3, SNRDE, mom, 
error) ; 


plot (SNRDB,error, 'y+') 

hold on 

end 

fclose(fpl); 

%* To calculate the dynamic range 
if rem(fsl1,2)==0 


x=[fs1/2 + fs2*fs3; fsi*fs2/2 
elseif rem(fs2,2)== 


oe 


To check whether fsl is even 

fs3;>"fsi*<ts3 + £S2) ; 

To check whether fs2 is even 
x=[f£s2/2 + fsi*fs3; fsl*fs272 fs3; fs2*fs3 + fs1]; 

elseif rem(fs3,2)== To check whether fs3 is even 
x=[fs3/2 + fs2*fs1; fs3*fs2/2 + fsi1; fsl*fs3 + fs2]; 

fsi,fs2 and fs3 are odd 

x=1/2* [fsl + fs2*fs3; fs2 + fsi*fs3; fs3 + fs2*fsl]; 


oe + 


+ 


oP 


oP? 


else 


end 
DR=min (x); 
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This function solves the general linear diophantine equation 


% m2*b1l - ml*b2 = 


k and returns the value bl 


function a=glde (m1,m2,k) 


° 


%* Calls function "lde" to calculate bl, b2 and na 


[b1,b2,na] =lde(m1,m2) ; 


oe 


ov? 


mult=k/na; 
if (k-mult*na) == 


bil=b1*mult; 
b2=b2*mult; 


mtest=bl1; 
mdl=m1/na; 
md2=m2/na; 
mx=b1; 
mx=mx+md2 ; 


To check whether the equation is solvable. 
na must be a factor of k for the equation to be solvable. 


%* Equation is solvable 


° 


* These new values solve the diophantine equation 


° 


ove 


To check whether bl and b2 are the least values 
that satisfies the diophantine equation 


oe 


while (abs (mx) -abs (b1))<0 


bl=mx; 


b2=b2-mdl1; 
mx=mx+md2 ; 


end 


if (mtest-bl)== 


mx=b1; 


mx=mx-md2 ; 
while (abs (mx) -abs (bl) ) <0 


end 
end 
end 


bl=mx; 
b2=b2+md1; 
mx=mx-md2 ; 


56 


This function solves the linear diophantine equation 

ml*b1l + m2*b2 = na where ml and m2 are the sampling frequencies 
and na is the greatest common divisor 

and returns the value bl, b2 and na 


oP oP oP o1? 


oe 


ml and m2 are assumed positive 


oe 


function [b1,b2,na]=lde(ml1,m2) 
% Initialize bol, bo2, bl and b2 


@1!= li, 
woz —U - 
ek—0 > 
lo eaL 


* Place ml and m2 in ma(dividend) and na (divisor) respectively 


ma=m1 > 
na=m2; 


% Calculate quotient and remainder 


iquot=fix(ma/na) ; 
irem=ma-na*iquot; 


% J£f remainder is not zero, reset dividend and divisor 


while irem>0 
bo3=bol1-iquot*b1; 
bo4=bo2-iquot*b2; 
bot=bar- 
bOZ-—D2. 
bil Dos ; 
b2=bo4; 
ma=na; * redefine dividend and divisor 
na=irem; 
iquot=fix(ma/na); % reapply Euclidean algorithm 
irem=ma-na*iquot; 


calculate new coefficients of ml and m2 


ov 


redefine bol, bo2, bl and b2 


oe 


end 


Si 
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Kex Ess TMS320C54x DSKplus 


Algebraic Assembler 








The 'C54x DSKplus is a low-cost design tool that gives designers a 


Code Explorer Debugger working knowledge of DSP code. From this foundation, designers can 
begin building complete 'C54x DSP-based systems. Priced at US $149, 

System Requirements the 'C54x DSKplus (part no. TMDS32000L0) is available from TI 
authorized distributors . 


How to Install 
The ‘C54x DSKplus builds on TI's industry-leading line of low cost, 


Bevond the DSKplus easy-to-use DSP Starter Kit (DSK) development boards. The 
high-performance board features the TMS320C542 16-bit fixed-point 
'C54x Software Support Files DSP. Capable of performing 40 million instructions per second 


(MIPS), the 'C542 makes the 'C54x DSKplus the most powerful DSK 
development board on the market. 


Other TMS320 DSKs include the 'C2x DSK . the 'C5x_ DSK , and'the 
floating-point 'C3x DSK . 











Key Features 


The 'C54x DSKplus includes: 


40 MIPS TMS320C542-based board 

TLC320ACO1 Analog Interface Circuit (AIC) 

'C54x DSKplus assembler, loader, Code Explorer debugger, and sample programs (3.5" disks) 
TMS320C54x CPU and Peripherals Reference Guide 
T™MS320C54x Algebraic Assembler Instruction Set 
TMS320C54x Datasheet 

TMS320C54x DSKplus User's Guide 

TLC320ACO01 Datasheet 

PC connector cable and universal power supply included 

US $49 discount coupon toward the purchase of the ‘C54x EVM 





-.°) 9/8/97 5:10 PM 




















DSKplus Key 
RaxtPaces : Benefits : 
TMS320C542 DSP (40 High-performance, very efficient architecture requires fewer 

MIPS, 16-bit) ‘| MIPS than competing DSPs to implement most algorithms. 
An easy-to-use, true Windows-based interface. Supports 
| Code Explorer debugger  { symbolic debugging, breakpoints, graphical animation, variable 
| interface | watch windows, file I/O, algebraic/mnemonic disassembly, 
| ‘| On-line help. 
| Symbolic debugging (Code Enables easy programmability by using labels for referencing 

Explorer) i} constants, variables, matrices by name. 


ee ae 


etecare 
Sweet at tease 6 sae 


Bypasses learning new DSP mnemonic instruction set specifics. 
Makes coding easier and more straight-forward. Easy one-step 
assembly and linking process. 


Algebraic assembler 


eh By 8h gy ORG wee se tate 
Sl Ge ee 


Demo programs / 


l : 
Application code Helps users get up-to-speed quickly : 


| 
' 3 
a. 
| 
| 
| 
i 
| 


Pateare s teases Get eee 
ee 


TLC320A.C01 Analog ‘| Low power dissipation, 14-bit linear resolution, programmable ! 


: sampling rates, anti-aliasing filter, and input gain; selectable 
| peace Chip auxiliary input; data read-back 
| 


« 
2 
a 








4 Socketed Programmable : Allows experienced designers to reprogram the PAL and change 
1 Array Logic (PAL) ii the way the host port interface works on the C54x DSKplus. a 
| 
| Universal Wy & Allows for immediate use out of the box; ideal for powering 
em iccaced ca -| daughter cards; filtered and regulated - thus no need for 


on-board voltage regulation. ; 





'C54x Algebraic Assembler 


The C54x DSKplus includes the algebraic assembler that speeds the initial code development process. The 
algebraic assembler does not require new users to learn a new DSP mnemonic instruction set, making 
coding easier and more direct. The assembler also utilizes a one-step assembly and linking process to 
simplify code debugging. The-software accomplishes this by using special directives to assemble code at an 
absolute address. 


Some extremely useful features include: 


® In-line Assembly expression analysis allows the assembler to work when defining complex variables 
or bit locations. 
Symbolic Debugging allows the user to reference variables by name instead of the physical address. 
Assembling conditional blocks of assembly code using .if/.else if/.end if directives. This is especially 
helpful when you want to conditionally assemble code via a command-line argument of internal 
assembly variable. 

© Support of .sect, .bss, .usect, .text, and .data sections. 





Code Explorer Debugger 


The 'C54x DSKplus debugger was developed by GO DSP Corporation in an effort to provide the first true 
Windows-based debugger for a DSK. The Code Explorer debugger supports debugging, a new feature 
available only on the DSKplus that allows the user to specify labels for referencing constants, variables, 
and marticies by name. Also, the debugger desktop environment is fully configurable and loaded upon 
entry into the debugger. This means that optional colors, fonts, and window sizes can be changed within 
the debugger and saved upon exiting. 


Some additional features of the debugger include capability of connecting files as /O, graphical animation, 
and data memory viewing. The file I/O capability enables users to connect files as inputs or outputs to any 
location within your application code. Therefore you can simulate different input sequences and data 
streams without having to physically generate them. 


Graphical animation allows you to view data in a graphical format, either with time domain or frequency 
domain and in a variety of variable sizes (i.e. 8-bit signed char, 8-bit unsigned char, 16-bit, 32-bit, etc). 


Disassembly Window 
















The disassembly window displays the DSP code in : 





s.Dis-Assembly 





algebraic instructions. The variable names and pean ote A t+ e7ih 
alee — oa oo cae lal 8 iF ee wiles: = 
Si address is the first column andthe : stort ae 
aoe in column 2, :2229 F974 cal) SCOIINIT a 
machine code for the instructions are in column 2. ‘9¢92 9731p ner(PHST) = Sledh = 
indi location where the :6504 77138 mur(SP) = #0ftah ee 
Tate alleys ide ote os 0506 7712 mar(AR2) = #1200h = 
DSP program counter (PC) points. -0508 7092 SA MEME ole cseontsy Site Crtane 
pe aay = 
° . ° Cc = 0 ~ eS 
The disassembly window properties can be '059D wait 
050D Fa9s nop 
accessed by placing the cursor in the disassembly ‘OS0E F073 ew hes 
window and right-clicking and then choosing ‘0510 XINT ae 
properties. The disassembly window can display =j2i} £230 ae mer (te 
code in algebraic or mnemonic formats with direct {0513 7092 SAR29 © 25 
and are eee values shown in hex, 0816 6012 nar (31B) eee Res) heer 
decimal and even binary. :0S1 0 it (TC) goto restrt 
:0$14 FSEB return_enable 
EIB 2212" mar(AR2) = #1200h 
Data Memory Window lal seers ret thegtee Pes 


Fae ay Seve, a ane 
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The data memory window can be modified or 
replicated as needed. By placing the cursor inside 
the data memory window and right-clicking and 





aM parainw Options. - . 


“+ Data ical “Tile = enon 







then choosing properties, the user can change the | 0x0800: 
title of the window, starting address and even data |: Ay 
organization in the window. Valid display formats | 0z080C: 
include 8-bit signed/unsigned char, 00810: 
signed/unsigned long, floats, and others. The page |: 98814: 
ee ie 0x0818: 
field can specify either Data or Program memory | 0=081C: 3 
spaces. . 0x0820: Mp avert tse . Float ccs 
0x0824: = Se a ee “SESS 
'C54x CPU and Peripheral — 
Registers 
The two register windows in the 'C54x Code oe ina = 
Explorer debugger are the CPU and Peripheral : - 0500 GS? ~ OL?L <7 eT | 
. = O00:00G ARO = 10FA oat 
Registers. The C34x CPU Registers is the B - FEFEFFFFFD BP] - 1002 OVA - 0 
collection of registers which control the operation . T = FrFoD AP? = 1005 oOYB=0 
; AR2 ~ EFEF OVH - 1 
of the DSP CPU. The program counter, status ree ae ea 
register, and configuration registers are contained . STO = 1860 APS ~ o000 Cib «0 
within this window. Notice that bit values within ae a Aare = = 3H FROT = 0 55 
. =O ts r ¢ 4PT wv 
the register are brought out separately to make + DP = 9000 BX = FFEF CPL = 
modification and monitoring easier. ‘ ASH = 000C ARS = ARO = = : 
: : ; BRAF = 0 INTSM = (C=O 
The second window is the Peripherals window. bo = FEFF TMR = & i 


This window includes the registers for configuring it * 
the DSP peripherals like the serial ports and 





- : : : : , 2 = } 
timers. Modifications to this register can be done sae : Go hen 
by clicking on the register in the Peripheral i TCR = cocd spe 

[TRC¥ - 0000 TDER 


Registers window. 
Graphical Windows 


Graphical windows are extremely useful | ESF 
when trying to view a value of a register, (22525 
variable, or buffer. The graphic window [Fit -22| 
allows the user to animinate any value in [esyaeing 3 ) 
either data or program DSP memory. iiraconee: 
This is accomplished by placing a 
breakpoint anywhere in the application 
code and pressing the Animation button. 
Each time the DSP reaches the 
breakpoint the graphical windows are Rep yt eee Bae sca peekt Beth aed 
updated and refreshed. Haiaen ons heres HIP Bua ie ss 


— = a tHe TS 
Se pce 





The options window contains the 
graphics setup for the window. For 
example, the title can be changed to 
reflect the data being animated, the 
display buffer length can be changed, or 
the data read from the DSP can either be EEE - oa = 

a single value from a list (buffer) of values in Bree ee Or program memory. ie the e sampling Tate can 
be modified for correct displaying of the frequency data (FFT). The display can be viewed using 8-bit 
signed/unsigned chars, ints, long, floats, and even a log can be performed on the displayed data. 


4 
vom. 


Setting Breakpoints 
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| 


A breakpoint can be selected by either 
double clicking on a line in the 
Disassembly window or by Selecting the 
DEBUG-BREAKPOINTS in the Pull 
Down Menu. The Pull Down Menu will 
prompt you with a menu listing all the 
available symbols in the Symbols box. 
You can either select a breakpoint from 
the list of Symbols or by entering an 
address in the Address field. 


J 
\ 
eee er 


The Breakpoint dialog box contains the 
following fields: Address, Symbols and 
Breakpoints. If the address of the desired 
breakpoint is known, simply enter the 
value in the Address field. The Symbol 
field contains the list of all the symbols in the program. If the location address of the breakpoint is labeled, 
simply type the label name and press add. 
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Setting Probe Points 

Probe points allow the update of a 2% , 

particular window or the reading/wniting Se 

of samples from a file to occur at a N eaceee ae e ee 

specific point in an algorithm. This kees Dees 

effectively "connects a signal probe" to iE ee  . : 
that point in the algorithm. Ses B 
When a graph window object is created, 4 -pabsteeie ers BA Re ee 


it assumes that it is to be updated at every if: 
breakpoint. However, this attribute can be jf 
changed and the window can be updated | 
only when the program reaches the 
connected probe point. After the probe 
point is hit, and the window is updated, 
execution of the program is continued. 
This optimizes the display of the graph 
window and also allows you to keep a history of the signal even when the data on the DSP is not valid. 


i 

H 

+? 
iF 
tp: 
IE 
iF 


EE Ts FD RT eee ey oe, 


0s ge em sess 


With the combination of Code Explorer's File I/O capabilities, probe points can also used to connect 
streams of data to a particular point in the DSP Code. When the probe point is reached in the algorithm, 
data is streamed from a specific memory area to file, or from the file to memory. 


Using File I/O 
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Code Explorer allows the user to stream 
data onto (or from) the target from a PC 
file. This allows the user to simulate code |}'= 
using known sample values. Note that 4 
this file I/O feature is not intended to 
satisfy real-time constraints. The File 
Input/Output feature uses probe points. 
When the execution of the program 
reaches a probe point, the connected 
object, whether it is a file, graph or 
memory window, is updated. Once the 
connected object is updated, execution 
continues. Using this concept, if a probe 
point is set at a specific point in the code 
and then connected to a file, file I/O functionalities can be implemented. 
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System Requirements 


A 386, 486, or Pentium PC with a 3.5" disk drive 
4-bit parallel and/or 8-bit bidirectional parallel ports. 
A minimum of 4Mbytes of memory 

Color VGA monitor 

Windows 3.1 or Windows 95 

ASCU editor 





How to Install 


When connecting the DSKplus to your PC, it is highly recommended you turn 
off your PC's power to make the connections below: 


Connect the DB25 cable (female) to the PC's Parallel port (male). 

. Connect the DB25 cable (male) to the DSKplus board (female). 
Connect the power cord (NEMA cable) to the 5 volt power supply. 
Connect the 5-pin DIN-to-5.5mm adapter to the power supply’s 5-pin 
DIN connector. 

Plug the power supply power cord to the wall outlet. 

Plug the 5.5mm connector into the power jack of the DSKplus board. 





PY oo 


oe 


At this point the green power LED is illuminated and power is supplied to the 'C54x DSKplus board. If the 
Green LED is not illuminated, check the connections on the power supply and power cord. 


Installing the software 
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The DSKplus kit includes two 3.5” floppies labeled Disk #1 and Disk #2. To 
install the software correctly, please follow the steps below: 


I. 
2. 


3h, 


Insert Disk #1 into the 3.5" drive. 

From the start menu (Windows95) or the Files menu (Windows 3.1) 
select the Run.. option. Type a:\setup.exe 

The installation script will appear. You will be asked to select a 
destination directory. By default it will select the DSKplus directory. 
Enter the directory name if you would like to specify a different 
directory. 





. When prompted, insert Disk #2 into the 3.5" floppy drive. 
. When installation has completed, the installation will inform you that the installation was 


successful. At this point a Code Explorer Group will appear. 


Starting the Debugger 


To start the debugger, click on the tcon located in the Code Explorer Group or desktop. The Code Explorer 
background and windows will appear with the Setup Box shown active. 


Select the port which is connected to the DSKplus board. Ifforsome reason (QEyGRRRI 
the port is not listed, the port address can be modified by typing in the address Fet== 


int the text box. 


As a result of selecting the correct port and proper hardware connections, the =P Ackbess > 


debugger will fill its windows with data and the DSKplus 1s now functioning. 
If for some reason the debugger responds with the error "Can't initialize 
Target DSP", follow the directions in the error box. 





Troubleshooting 


1. 


Zs 


Is the power on? Be sure green LED is illuminated. If not, a loose power cable is hampering your 


setup. 
Is the parallel port cable connection secure? In many new DSKplus boards and parallel port cables, 
substantial pressure many be needed to connect the cables. Connect the cable to the DSKplus board 
by placing the thumb behind the DB-25 connector. Take the cable connector chassis and place 
between the index and middle fingers. Align the connectors and press the fingers together. 


. The port selected is not being "Captured" by Windows 95. Capturing is used by Windows 95 to 


allow DOS programs access to printers. The port can be released by going into the control panel 
and selecting the printers icon. Highlight any printer and go to the File pulldown on the command 
bar. Select properties and then the Details tab. The Details tab includes a button named End 
Capture... Click on this button and select the LPT port where the DSKplus board is connected. If 
the LPT is not listed, then the port is not captured (select cancel) and proceed to number 4. 


. The port selected is configured as an EPP or ECP port. The DSKplus board supports 4-bit 


unidirectional and 8-bit bidirectional parallel ports. The DSKplus does not support EPP and ECP 
ports. To check the port configuration, exit out and reboot your system. At the point where the 
BIOS Setup routine can be selected, press the keyboard sequence to enter the BIOS (usually 
CTRL+ALT+ESC). Confirm that the paraliel port is setup as '8-bit', ‘bidirectional’ or ‘standard.’ 
Specifically, not an EPP or ECP port. If problems persist, run the included selftest program. 








er ere ee ee) - .« ‘m+ 


Beyond the 'C54x DSKplus 


With highet performance than any other DSK available today, the 'C54x DSKplus offers a rich 
development environment for benchmarking and evaluating code in real-time. The 'C54x DSKplus is 
designed as an easy-to-use entry into the world of high-performance fixed-point DSPs. 
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However, as your design experience grows, you may require additional functionality and expanded 
capabilities. To meet these needs, TI offers a comprehensive line of evelopment tools for the TMS320 DSPs 
that support the design process from system concept to production. 


Other 'C54x Development Tools 
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APPENDIX C 


C AND ASSEMBLY LANGUAGE CODE FOR DSP HARDWARE 


—eKKKKRKKAKRAE EE EA RARER E AAA YX aX RRA ES 


; File: FirstAp1.ASM 
; a. buffer size here 
; b. A and B registers in 


d 


When sampling frequency is changed, need to change 


ACOlinil asm 


c. sampling frequency in dftsort.c 
d. buffer size in hostappl.cpp 


se KKEKKKKKKKEKEKKKKEKKKEKKEKKEEKEKEKKKKEKKKKEKEKKEERKEKEKEKEEEKEKKKEKKKEEKKEKKKKE 


.width 80 
.length 55 


.title "FirstApp program" 


.mmregs 


.setsect ".text", 0x500,0 
.setsect "vectors", 0x180,0 


ee 
— ee ee ee eee ee ee ee ee ee ee ee eee ee ee eee ee ee ee eee ee ee ee ee ee ee ee ee ee ce ce ee ee ee ee ee ee ee ee 


ce ee ee ee ee ee ee we ee Sw 


se ee ee ee 


cc cr cr cc cc ww we ew ee ee es ee ie oe ee 
ee Oe 


SISO Gubideic coisas 


SO ee _____. _._._ 


.copy "c:\dskplus\inits\vectors.asm" 


text 

Start: 
call ACO1LINIT 
pmst = #01la0h 


=e 


sp = #0ffah : 
ar2 = #1200h ; 
*ar2+ = data(#0bh) r 
imr = #280h 

intm = © : 


wait nop 
goto wait 


$000 te eee eee Receive 
AINT: 
b = trcev ; 
b = #0FFFCh & b 
*ar2+ = data(#0bh) ; 
ean =. 5 
TC = (@ar2 == #1471h) 


=e me ™e 


LE MTC) Goto restnre 
return_enable 
Bese G 


set up iptr 

init stack pointer. 

pointer to receive buffer at 1200h. 
store to rev buffer 


ready to rcv int's 


Interrupt Routine ------------------ 
load acc b with input 

store to rev buffer 

transmit the data. 


change here if fs changes 
stop if rev buffer is at 1471h 


69 


ar2 = #1200h ; §SCGemem bit ...no int's 
hpic = #0ah ; flag host task completed 


return_enable 


.copy "c:\dskplus\firstapp\acOlinil.asm" 

.end 
g EK KKKKEKEKKEKEKEKEKEKRKEEKRREREKRKREEKKEKREKEEEKEKEKEKEEKEKEEKEEKEKEEKKEKEKEKKKKKKKKKKKKK 
; File: FirstApp.ASM -> First Application program for the 'C54x DSKplus 
buffer size here 
A and B registers in ACOlini2.asm 
sampling frequency in dftsort.c 
. buffer size in hostapp2.cpp 


ano 


eK KKKKKEKKKKKKKKKKEKKK KKK KKKEKKKKEKKEKKEKKEKKKEKKEKKKKEKKKKKKEEKKEKKKKKKKKKKKKKKKK 


.width 80 

-length 55 

.title "FirstApp program" 
.mmregs 

.setsect ".text", 0x500,0 


.setsect "vectors", 0x180,0 


mc ew ec cc we ce ce rr cc cc cc ce ew ww we we we = 
— ep cc ce ae cet ccc cr crc cre i er rc cc cc ce cc ce ce ce ce ee ce ce ee ee ee ee ee ee ee ee ee 


Cn ee ee 
ee eT LL 


.sect "vectors" 
.copy "c:\dskplus\inits\vectors.asm" 


ZCSCXE 
State: 
call ACOLINIT 
pmst = #01a0h 7. Set Up Apter 
sp = #0ffah ; init stack pointer. 
ar2 = #1200h ; pointer to receive buffer at 1200h. 
*ar2+ = data (#0bh) ; store to rcv buffer 
imr = #280h 
imem, — 0 ; ready to rcv int's 


wait nop 
goto Wait . 


BSS SSS SSeS SSS Se Receive Interrupt Routine ------------------ 


AINT: 
meine 7 ; load acc b with input 
b = #0FFFCh & b 
*ar2+ = data (#0bh) peStore werrev buifer 
[gob.gg = Je) ; transmit the data. 
TC = (@ar2 == #01400h) ; change here if £s change 
if (TC) goto restrt ; stop if rev buffer is at 1400h 


return_enable 


70 


EEStre 


ar2 = #1200h ; Seteintm bit ...nO ime 
hpic = #0ah ; flag host task completed 
return_enable 

$0 ere terre ------- end ISR ------------ rrr creer creer 


.copy "c:\dskplus\firstapp\acOlini2.asm" 
.end 


71 


a oe oe oe oe oe ee oe oe a a ee ee ee a a a 2 a a ae 
f a 


; 


: File: 


« 
é 


ACO1LINI1.ASM -> ACOl Initialization Routine 


aK KEKE KEKEKEKEKEKEKKEKEEKKEEKEEKEEKEKEEKEKE KEKE KKEKKEKEKEKEKKEKEKEKEKKKEKKKEKEKEKEEKEKEKKEKKEKKEKKK 
‘ 


* + + + + 4 + + © © £ &£ & £ FF FF H FF FF FF HE HF FF FF F F HF HF HK F HF 


.width 80 

-lenoth 35 

.title "ACOl Initialization Program" 
.mmregs 


REGISTER 


0000 
0000 
0000 
0000 


0000 


0000 


0000 
0000 


(binary) 


0000 
0000 
0000 
0000 


0000 


0000 


0000 
0000 


REGISTER .set 
REG1 SCE 
REG2 .set 
REG3 cee 
REG4 .set 
REGS -set 
REG6 -set 


0000 
0000 
0000 
0000 


0001 


0010 


0100 
1000 


Any combination of 
number to the REGISTER constant. For example to initalize Registers 4 
and 5, REGISTER = 
initialization is included in the ACO1INIT module. 
module will load REG4 and REGS values into internal ACOl registers. 


Obh 
1feh 
Zale’ 
300h 
40dh 
50 tr 
600h 


— 
= 


0001 
0010 
0100 
1000 


0000 


0000 


0000 
0000 


The constant REGISTER should be set to 
register: 


initialize 
initialize 
initialize 
initialize 


=-> 
SS 
-> 
- > 
initialize 
initialize 
initialize 
IM alls Ze 


=— > 


wom 


KEKKKKKEKEKKEKEKEKEKEKEKKEKEKEKEKEKEKEKEKEKEKEEKEKEKEKEEKEKKEEEKEKEKEKEKEKEKEKEKEKEKKEKEKKEEKEKEKEKEKKEKEKKKKEKEKEKE 


Certain ACO1l registers can be initialized using a conditional assembly 
constant. By setting the constant REGISTER to the appropriate value, 

the assembler will either include initialization for certain registers 
or ignore register initialization. 


include the following ACO1 


Register 
Register 
Register 
Register 


Register 
Register 


Register 
Register 


registers can be initialized 


PWN eB 


7 
8 


(A Register) 

(B Register) 

(A' Register) 
(Amplifier Gain- 
Select) 

(Analog 
Configuration) 
(Digital 
Configuration) 
(Frame-Sync Delay) 
(Fram-Synce number) 


by adding the binary 


18h. Upon assembly, only code for register 4 & 5 


° 
é 


Powerup default values: 


ez 


Haley! 
212i 
300h 
405h 
S0in 
600h 


When called the 


Register 4 is always loaded to get a 6db input gain. This sets full- 
scale to 3v(p-p input) due to the single-ended AC0O1 configuration. 


REG7 .set 700h ; 

REG8 -set 80ih : 

ACOLINIT: 
ema 0, ; reset 
Behe = AL ; 
(ere FIL lOle) j 
imr = #280h ; 
tspc = #0008h : 
taxos—- 0h ; 
tspc = #00c8h ; 
Soe IL ; 


oem eee eee ee 


; Register init's 


.eval REGISTER & 1h, SELECT : 
fib SeEBel fae ny ; 
a = #REGI1 : 
call REQ2 : 
.endif 

-.eval REGISTER & 2h, SELECT : 
-if SELECT = 2h 

a = #REG2 

call REQ2 

.endif 

.eval REGISTER & 4h, SELECT , 
Pee Sh bber = 4h 

a = #REG3 

call REQ2 

.endif 

-eval REGISTER & 8h, SELECT . 
-if (SEEECT — son 

a = #REG4 

call REQ2 

.endif 

.eval REGISTER & 10h, SELECT =; 
-if SELECT = 10h 

a = #REG5 

call REQ2 

.endif 

-eval REGISTER & 20h, SELECT ; 
,TESSSEECT =-20h 

a = #REG6 

call REQ2 

.endif 
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700h 
80ih 


ac0Ol 


stop timer 
wakeup from idle when TDM Xmt int 
stop TDM serial port 

send 0 as first xmit word 
reset and start TDM serial port 
release acOl from reset 


Gisable all int service routines 


if REG1 then include 


load Acc A with REG1 
Call REQ2 subroutine 


if REG2 


if REG3 


if REG4 


if REGS 


if REG6 


then 


then 


then 


then 


then 


include 


include 


include 


include 


include 


this 


value 


this 


this 


this 


this 


thas 


source 


source 


Source 


source’ 


Source 


SsOuUrRCEeE 


REQ2 


.eval REGISTER & 40h, SELECT ; if REG7 then include this source 


,if SELECT 
a = #REG7 
call REQ2 
.endif 


.eval REGISTER & 80h, SELECT ; if REG8 then include this source 


,2tE SELECE 
a = #REG8 
call REQ2 
.endif 
return 


ifr = #080h 
tdxr = #03h 


idle (1) 
ceaxr = 4 
ifc = #080h 


idle (1) 
oboe == 16)al 
mer = +08 0h 
idle (1) 
return 

.end 


clear flag from IFR 
request secondary when ACOl starts 


wait for primary to xmit 
send register value to serial port 
clear flag from IFR 


wait for secondary to xmit 

send neutral state in case last init 
clear flag from IFR 

walt for neutral state to xmit 
return from subroutine 
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—KEKKKKKKHKEKKEKEKKKEKKKKKKKKEKKKEKEEKEKREREEEEKEKEEKEKKKKEKEEKEKEKKKEKKEKEEEKKEKKEK 
P) 


2 
s 
e 
s 


c 


File: ACO1LINI2.ASM -> ACO1 Initialization Routine 


KKK KKKKKKKEKEKKKEKEEKKKEKKERKKEKEEKEKKEEKEKEKEKKEEKEEEEKEKKEEKKEKKEKKKKRKRKKRKKKESE 
, 


+ %* + + + + + £ FF €* F&F € F FF F HH FF F&F FF FF FF HF HH FF FF FF F HH HF HF FF HF 


Widin 80 

J lenges 55 ; 
title “ACO] Initialization, Program" 

.mmregs 


KKK KKKKEKKKKKEEEEEEKEKKKKEEKEKEKEKKEKEKKEKEKKEEEKEKEKEKEKEKEKEKKKKEKEKKEKEKKEKKKKEK 


Certain ACOl registers can be initialized using a conditional assembly 
constant. By setting the constant REGISTER to the appropriate value, 
the assembler will either include initialization for certain registers 
or ignore register initialization. 


The constant REGISTER should be set to include the following ACOl1 
register: 


REGISTER (binary) = 


0000 0000 0000 0001 -> initialize Register 1 (A Register) 

0000 0000 0000 0010 -> initialize Register 2 (B Register) 

0000 0000 0000 0100 -> initialize Register 3 (A' Register) 

0000 0000 0000 1000 -> initialize Register 4 (Amplifier Gain- 
Select) 

0000 0000 0001 0000 -> initialize Register 5 (Analog 
Configuration) 

0000 0000 0010 0000 -> initialize Register 6 (Digital 
Configuration) 


0000 0000 0100 0000 -> initialize Register 7 (Frame-Sync Delay) 
0000 0000 1000 0000 -> initialize Register 8 (Fram-Sync number) 


Any combination of registers can be initialized by adding the binary 
number to the REGISTER constant. For example to initalize Registers 4 
and 5, REGISTER = 18h. Upon assembly, only code for register 4 & 5 
initialization is included in the ACOLINIT module. When called the 
module will load REG4 and REGS values into internal ACO1 registers. 


Register 4 is always loaded to get a 6db input gain. This sets full- 
scale to 3v(p-p input) due to the single-ended ACO1 configuration. 


REGISTER .set Obh ; Powerup default values: 
REG1 -set 1feh * 112h 
REG2 -set 23Cn al 212h 
REG3 -set 300h ; 300h 
REG4 -set 40dh ;* 405h 
REGS -set 5 Odin ; 50ih 
REG6 -set 600h ; 600h 


@S 


REG7 .set 700h ; MOON : 

REG8 see 801ih : 80ih 

ACOIINIT: 
pee = ; reset acOl 
a evey se ll ; disable all int service routines 
ECE = Fon ; stop timer 
imr = #280h ; wakeup from idle when TDM Xmt int 
tspc = #0008h ; stop TDM serial port 
tdxr = #0h ; send 0 as first xmit word 
tspc = #00c8h ; reset and start TDM serial port 
soe et il ; release ac0Ol from reset 


.eval .REGISTER & 1h, SELECT 
-1£ SELECT = lh 

a = #REGI1 

call REQ2 

.endif 


.eval REGISTER & 2h, SELECT 
Sif SOLRECT =e2n 

a = #REG2 

call REQ2 


endif 


.eval REGISTER & 4h, SELECT 
.1f£ SELECT = 4h 

a = #REG3 

call REQ2 

-endif 


.eval REGISTER & 8h, SELECT 
Mae SeLECT = 8h 

a = #REG4 

call REQ2 

.endif 


.eval REGISTER & 10h, SELECT 
Sli Sseueec) = 10h 

a = #REGS 

call REQ2 

.endif 


.eval REGISTER & 20h, SELECT 
.1f SELECT = 20h 
a = #REG6 


° 
f 
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— ee ee ee eee eee ee i eee i 


if REG1 then include 


load Acc A with REG1 
Call REQ2 subroutine 


slat 


ahie 


ade 


acc 


nLae 


REG2 then 


REG3 then 


REG4 then 


REGS then 


REG6 then 


include 


include 


include 


include 


include 


this source 


value 


this source 


this source 


this source 


this source 


this source 


REQ2 


call REQ2 


endif 

.eval REGISTER & 40h, SELECT ; if REG7 then include this source 
.i1f SELECT = 40h 

a = #REG7 

call REQ2 

.endif 

.eval REGISTER & 80h, SELECT ; if REG8 then include this source 
Rane. SPR SGik ses /sl0)o 

a = #REG8 

call REQ2 

PenGanr 

ie elt ied) 

ifr = #080h ; clear flag from IFR 

bax — +03h ; request secondary when ACOl1 starts 
idle (1) ; wait for primary to xmit 

Eds = 7a ; send register value to serial port 
ifr = #080h ; clear flag from IFR 

idle (1) ; wait for secondary to xmit 

Saxe = 20h ; send neutral state in case last init 
Pir — -080n ; clear flag from IFR 

idle (1) ; wait for neutral state to xmit 
return ; return from subroutine 

.end 


T@ 


iment EER aie a RRS EAR RRR ce OS BO ee 
/* 

/* File: HOSTAPP1.CPP Source code for host application 

/* 

[RAK KRKKKKKKKK EEK KEKE KKK KKK KK KKK KEKE KEKKEK KEK EKEEEEKEKKERKEKRKEKREKEKKEKKKK EK | 
#include <H1I54X.H> 

#include <stdio.h> 

#includes<cOnloen= 

#include <stdlib.h> 

extern int daeameg|],statreg(],ctrilreg([]; 

extern int pport, portmode,Readdelay; 


void main(void) 


{ 


FILE *fp; 
if ((fp=fopen("datal.dat", "w") ) ==NULL) /* Open file */ 
{ 
‘elcecr (); 
printf("Cannot open file .\n"); 
exit (0)¢ 
} 
portmode=0; /* 4-bit mode */ 
Readdelay = 20; /* In case host slow*/ 
Gilteseia se; /* Clear the screen */ 


if ((pport=locate port()) >= 5){ /* Find the port. */ 
Denet "NO connection \n"); /* If no connection */ 


backout () ; /* then leave board */ 
exit (0) ;} /* in known state */ 
else{} 

_setcursortype (_NOCURSOR) ; /* Hide text cursor */ 
set _latch(1,1); /* Keep DSP running */ 
imt word =0, col=0; /* and bring PAL out*/ 
/* out of Tri-state */ 

col=0; 
qecess (1,1); /* go to home Lh // 
send_word(0x0808, C_SEND) ; /* Clear the HINT * / 
HEINE (10000) - /* Wait for nxt HINT*/ 
send_word(0x1200, A SEND) ; /* Goto 0x46 entries*/ 


/* before buffer x / 
for(int buf=0 ; buf < 0x271; buf++) 
/* change here if fs change*/ 


{ 

word = read word(D READ); /* Read word from pp*/ 
PEinth ('s4.4= "| word); /* Preinteit to scr °*/ 
EPrinee( tone sd \n", word), /* ,Outpucrto file * / 
if (Cotes —.13) { /* in 14 columns x / 


col=0; 
jenenloghese (uae alu) 
else{col++; } 


_setcursortype (_NORMALCURSOR); /* Ret normal cursor*/ 
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fclose (fp) ; /* Close file * / 
backout (); 
/* Leave board in * / 
exit (0); /* known state x / 


ae) 


[RR KKK KKK KKK KKK KKK KK KEK KK KKK KKK KK KK KK KKK KKK KKK KKK KK KKK KKK KKK KKK KKK KKK KKK / 
/* 

/* File: HOSTAPP2.CPP Source code for host application 

/* 

[RRR KRE KKK KKK EK KERR KK EEK KKK KEKE ERK KKK KKK KK EKEKKE KKK EK KRKEKEKKKKEKEKEKEKKE / 
#include <HI54X.H> 

#include <stdio.h> 

#include <conio.h> 

#include <stdlib n> 


e 


extern int datareg[],statreg[{],ctriregI[]; 
extern int pport, portmode,Readdelay; 


void main(void) 
{ 
FILE *fp; 
i= ((Ep=fopen("data2.dat", "w") ) ==NULL) /* Open file */ 
{ | 
elrscr(); 
Beumtcr ("Cannotvooen tlle -.\n")’; 


exit (0); 
} 
portmode=0; /* 4-bit mode */ 
Readdelay = 20; /* In case host slow*/ 
clrser or /* Clear the screen */ 
if ((ppert=lecate pert()) >= 5){/* Find the port: */ 
PpeintmlueNe CONnnece1on\n")>; /* LL no connection 27 
backout () ; /* then leave board */ 
ext (0). | J* in known ceatbe * / 
else{} 
_setcursortype (_NOCURSOR) ; /* Hide text cursor */ 
set_latch(1,1); /* Keep DSP running */ 
int word =0, col=0; /* and bring PAL out*/ 
/* out of Tri-state */ 
cen); 
Setocy (4): /* go to home x [ 
send_word(0x0808, C_SEND) ; /*@@lear the HINT a 
HINT (10000) ; /* Wait. for asses arias 
send_word(0x1200, A_SEND) ; /* Goto 0x46 entries*/ 


/* before buffer x / 
for(int buf=0 ; buf < 0x200; buf++) 
/* Change here if fs is changed*/ 
{ 


word = read _word(D READ); /* Read word from pp*/ 


Pept r( 44.45  eword); /* Pointeit to scr */ 

Porm Eo od \n"; Word) >; 7*seuteput to file * / 

if (col >= 13) { /* in 14 columns * / 
eou- Ur 


Breroer <i \n!).; | 
else{col++; } 
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_Ssetcursortype (_NORMALCURSOR) ; /* Ret normal cursor*/ 


fclose (fp) ; /* Close file * / 
backout () ; /* Leave board in * / 
exit(0); /* known state * / 


Sa 


[8K KK KK KR KR RK KK KK KR RK KR RR RR RK KR RK KR RK KK KK KK / 


/* 
/* File: 
/* 


MAIN.C Source code for main program 


[RR KR RK KKK KKK KK KK KK KK KKK KKK IK KK KEI KKK KKK KKK KKK KK KKK KKK KKK KR KKK KK KK / 


#include 
#include 
#include 
#include 
#include 


#define 
#define 


main () 


{ 


i qiks weg ale 
"Stdde£ say 
"Stdliban’ 
"Stdto wi 
" COnd@ an — 


fsi 625 
fs2 324 


FILE *fpi, *fp2; 

ieee. G2, xll[fsil, xZlts2); 
Miles, OUL, Lreq; 

int gide (int k); 


Geulble.xrOollislijwx2ol (isil> xo1 (Esil- 
float .Oi= 932145926, tpi; 

Dnt ries 

int 25 ease 

double maxl1, max2; 


double xroZ2 [isZivex1oz (fs2], xo2lfs2!; 
int Jegeafke2. 2; 


if ((£pl=fopen("data5.dat", "rt") ) ==NULL) 
{ 
Gimacr (); 
pmimea.* Cannot Open file .\n") ; 
exit (0); 


} 


if ((fp2=fopen("dataé6é.dat", "rt") ) ==NULL) 
{ 
Girser()> 
printf ("Cantiot open file .\n"); 
exit (0) ; 
} 
for (n=0;n<fsl1;n++) 
{ 
ESGene(tpl, "se ", &di); 
>it (gil se ok 


} 


for (n=0 ;:n<fs2;n++) 


{ 


Bocaire p2, "se ", &d2) ; 


314 


/* Open file */ 


/* Open file */ 


x2 (nN) =G2; 


} 


fclose (ipl 
fclose (fp2z)- 


tpil=Z-ea 

for (u=0;u<fs1;u++) 

{ 
xEOl (ui —02 05 
xiol [u] =0.0; 
for (n=0;n<fsl1;n++) 


{ 

/*-- Xr[u] = (1/fs1) sum {xr[n].cos(2PI.u.n/fsi)} --*/ 
xrol{u] = xrolf{u] + x1[n]*cos(tpi*u*n/fs1) ; 

/*-- Xifu] = - (1/fs1) sum xr[n].sin(2PI.u.n/fs1) --*/ 
xiol[u]. = xiol[u)] = x@iai*sin(teri-n/ tomy 

} 


xrol {uj =xrol1 [fu] /fs1; 
xiol [u] =xiol [u] /fsl; 
sod fuil—eemet (xrol [ul *xrol (ul +xieistal sro inl 
} 
adfiti=0; 
nl=fs1/2+1; 
max1l=xo1{[0]; 
for (i=l;i<n1;i++) 
{ 
if (xol{i] > max1) 
{ 
aftl=i1; 
maxl=xol [i]; 


} 


for (u=0;u<fs2;u++) 
{ 
Meo (| — 020; 
Sez =O 0; 
for (n=0;n<fs2;n++) 
{ 
xro2[{u] = xro2[u] + x2{[{n]*cos(tpi*u*n/fs2) ; 
xio2[{u]) = xio2[u] = x2[n]*Sin(eps*oew csZ)- 


xro2 [u] =xro2 [u] /fs2; 
xio2 [u] =xio2 [u] /fs2; 
xo2 {u] =sqrt (xro2 [ul] *xro2 fu] +x102 [u]|/*xie27u)]) - 


} 


dadft2=0; 
n2=fs2/2; 
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glide (k) 


/* This section solves the linear diophantine equation fsl*bl + fs2*b2 


max2=xo2 [0]; 
for (j=1;j<n2;J3++) 


{ 
Lf (xe2zg) > max2) 
{ 
chege=] 
max2=xo2 [j] ; 
} 
} 
K=dEtZ-dEt1;> 
out=glde (k) ; 


freq=abs (dftli+out*fs1) ; 


printf("The frequency is ¢d", freq); 


far?) ; 


float mult,b1,b2,mtest,mdl,md2,mx; 


na where fsl and fs2 are the sampling frequencies 


and na is the greatest common divisor and returns the value bl, b2 and 


na. fsl and fs2 are assumed positive */ 


float bol,bo2,ma,na,irem,bo3,bo4; 
Ine / LQUCE - 


okie. 
boz— 
Or 
b2=1; 


/* Place fsl1 and fs2 in ma(dividend) and na 


ma=fsl1; 
na=fs2; 


/* Calculate quotient and remainder */ 


/* If remainder is not zero, 


iquot=ma/na; 
lrem=ma-na*iquot; 


while (irem>0) 


{ 


(divisor) respectively */ 


reset dividend and divisor */ 


bo3=bol1-iquot*bl1; /* calculate new coefficients */ 


bo4=bo2-iquot*b2; 


bol=bl1; /* redefine bol, bo2, bl and b2 */ 
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bo2=b2; 


eli —lele)s) 

b2=bo04; 

ma=na; /* redefine dividend and divisor */ 
na=irem; 

iquot=ma/na; /* reapply Euclidean algorithm */ 


i1rem=ma-na*iquot; 


/* To check whether the equation is solvable. na must be a factor of k 
for the equation to be solvable. */ 


mult=k/na; 
if ((k-mult*na) ==0) /* Equation is solvable */ 
{ 
bi=b1*mult; /* These new values solve the * / 
/* diophantine equation * / 
b2=b2*mult; 
mtest=bl1; /* To check whether bi and b2 * / 
/* are the least values that * / 
mdi=fsli/na; /* satisfies the diophantine equation */ 
mdad2=fs2/na; 
(eee 


mx=mx+md2; 


while ((abs(mx)-abs(b1)) < 0) 


{ 
bi=mx; 
b2=b2-mdl1; 
mx=mx+md2 ; 


} 


if ((mtest-b1) ==0) 


{ 
mx=b1; : 
mx=mx -md2 ; 
while ( (abs (mx) -abs (b1) ) <0) 


{ 
bi=mx; 
b2=b2+mdi1; 
mx=mx-md2 ; 
} 
} 
mecurn ( (int) bl) > 


} 
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